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Why Transactions? 


ERNST WEBER 
Microwave Research Institute, Brooklyn, N. Y. 


With the growth of professional organizations, 
the publications program becomes inevitably a 
serious problem. IRE has been no exception in this 
respect and publication policies have been the cen- 
ter of many subjective and objective discussions. 
Subjectively, each member would like to see the 
main organ of IRE carry important news and read- 
able articles covering fully his own particular field 
of interest. Objectively, the Editorial Department 
has endeavored to satisfy the necessarily widely 
divergent individual desires while also trying to 
maintain the very high standard of papers within 
rigidly confining space limits. It should be obvious 
that any ideal solution is unattainable through the 
PROCEEDINGS alone. 

With the initiation of the professional group sys- 
tem, a solution of the difficult publication problem 
appears to be at hand in a most natural manner. 
Each professional group may publish in its trans- 
actions almost at a moment’s notice any paper or 
group of papers of interest to its circles, without 
long formal reviews, painful surgeries, and un- 
happy feelings of authors. Moreover, transactions 
are destined to become the repository for outstand- 
ing papers in the particular field of each profes- 
sional group, contributing collected papers on 
specific subject matter! Surely, it is now the re- 
sponsibility of the professional groups to maintain 
the standard of the papers by effective use of their 
own editorial prerogatives. What was manifestly 
an impossible task for the PROCEEDINGS has be- 
come possible by delegation to 22 professional 
groups; a full, efficient, rapid, up-to-date reporting 
of scientific and technological developments to 


- each interested sector of the IRE membership! 


There might be a feeling on the part of some au- 
thors that the circulation of his contribution might 
be restricted by publication in the TRANSACTIONS 
as compared with publication in the PROCEEDINGS. 
This is not justified on at least two counts: unless 
the contribution is of such wide membership inter- 


| est that it warrants full distribution to all mem- 
| bers, it will not be published at all in the PROCEED- 


INGS; and the PROCEEDINGS will carry abstracts of 
all TRANSACTIONS papers and, in addition, will in- 
clude their authors and titles in the annual TRANS- 
ACTIONS index. Actually, therefore, the TRANSAC- 
TIONS will constitute specialized “sections” of the 
over-all Institute publications and thus rank with 
the PROCEEDINGS from the professional point of 
view. Indeed, the PROCEEDINGS will now be able 
to concentrate upon a single objective, namely to 
publish matters of wide interest, and preferably of 
interest to most members of IRE. 

Just what the division of subject matter between 
the PROCEEDINGS and between the TRANSACTIONS 
of the Professional Groups might be, whether, in- 
deed, a clear-cut division should be made, is yet 
an open question. Much will depend upon the use 
that Professional Groups will make of this un- 
paralleled opportunity for self-expression. One pos- 
sibility might be that the PROCEEDINGS include 
tutorial, informational, and basic scientific papers 
of high quality and of fundamental interest to a 
large majority of IRE members, granting the au- 
thors ample space for the full development of their 
ideas and thus making the papers valuable and 
lasting reference material. Most of the more spe- 
cialized technical papers would then go automati- 
cally to the TRANSACTIONS which might, of course, 
also carry news of plans and activities of the par- 
ticular Professional Group and its chapters. 

Surely a publication policy which can achieve 
satisfaction of the unifying membership interests 
in the PROCEEDINGS and of the divergent member- 
ship interests in the TRANSACTIONS of the Profes- 
sional. Groups is as close to an ideal solution of the 
intricate publication problem as we should hope 
for. Its success must, of course, depend upon the 
wholehearted support of each and every member 
of IRE and particularly upon the hard work of the 
laboring authors and of the editorial staffs. “By 
their deeds ye shall know them” . . . I feel sure the 
proper deeds will be forthcoming in the interest 
of the members, who after all constitute IRE, no 
more and no less. 
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news and views 


Tue Mucu ABUSED “DECIBEL” 


ISCUSSION revolving around the terms “decilog” 
1) and. “logit” has been going on for some time, and 

seems to be coming to a head. Despite the efforts 
of some stalwart purists, engineers have persistently 
used the “decibel” loosely, thereby causing painful er- 
rors. 

E. I. Green, Chairman of the American Standards 
Association subcommittee concerned with communica- 
tion terms, has proposed a new unit for logarithmic 
ratios to replace the overworked decibel. His thinking 
on this problem is incorporated in the ensuing article of 
his writing. 


A New UNIT For LOGARITHMIC RATIOS 


The use of a simple logarithmic unit based on power 
ratio for the distortionless measurement of telephone 
transmission was first proposed by R. V. L. Hartley. 
Such a unit, with the name of “transmission unit,” was 
adopted in 1924 by the Bell System in place of the previ- 
ous “ mile of standard cable.” Four years later the Inter- 
national Advisory Committee on Long Distance Teleph- 
ony in Europe adopted the “bel” and the “neper” as 
alternative units for transmission measurement and the 
transmission unit was rechristened the “decibel.” While 
the bel was the unit officially selected, the decibel has 
proved so much more convenient in size that the bel is 
now all but forgotten. The definition of the decibel is 
simple, the number of decibels corresponding to the ratio 
of any two powers P; and P» being equal to 10 logis 
Py Ps 

Through the years the decibel has found wide utility 
in communication and allied fields. The logarithmic ex- 
pression is handy because it saves explanatory words 
and because for large ratios it requires a minimum of 
figures. It is simpler to say that two powers differ by 50 
db than either to say that they bear a ratio of 100,000 
to 1, or to employ the exponential expression that one is 
10° times the other. More important than this, it is 
far easier to add and subtract db’s than to multiply or 
divide the corresponding ratios. 


So useful has the db proved that scientists have been 
lured into extending its meaning, in some instances with 
a degree of reason, but in others improperly, in ways not 
conforming to the original definition. It was rather natu- 
ral to extend the decibel to quantities related to power, 
for example, voltage squared or current squared. Thus 
the number of decibels corresponding to a scaler current 
ratio I,/I, or a scaler voltage ratio Vi/ V2 was taken as 
20 login Vi/ Ve, respectively. Initially the db was applied 
to voltage or current ratios only in cases where these 
ratios were the square roots of the corresponding power 
ratios, but it was gradually extended to voltage and cur- 
rent ratios where this was not the case, a practice which 
has received the sanction of defining bodies, provided 
specific statement of conditions is given. The decibel 
has also been extended, with approval of defining bodies, 
to ratios of sound intensity (using 10 Xlogio) and ratios 
of sound pressure and sound particle velocity (using 
20X logio) ; t 

While even this extension of usage tended to introduce 
a certain amount of confusion, the confusion has in re- 
cent years become more confounded as engineers have 
discovered that a logarithmic unit is convenient for ex- 
pressing ratios of quantities entirely unrelated to power, 
and in default of an accepted unit for this purpose have 
pirated the db. A number of people are beginning to 
use the db as a yardstick for ratios of values of length, 
or size, or frequency, or ohms, or permeability, or dol- 
lars, or whatnot. Usually they ignore the basic difficulty 
which arises from the fact that db’s are equal to 10 
times the logarithm of the ratio when applied to powers 
and 20 times the logarithm of the ratio when applied to 
voltages or current. Hence in using the db for other 
kinds of ratios, some people have made it always equal 
to 10 times the logarithm, some 20 times, and some have 
used either 10 or 20 times according to whatever similar- 
ity they found or fancied to exist between the quantity 
under consideration and either power or voltage. 

This extension of the db to magnitude ratios of any 
kind has been a matter of growing concern to defining 
bodies and others. Vigilance and deprecation, however, 
have only partially arrested the trend. To illustrate the 
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chaotic state of affairs which now exists, the writer, as 
Chairman of the ASA C42 Subcommittee on Definitions 
of Communication Terms, has been bombarded with 
letters divided about equally between pleas for exten- 
sion of the db to new kinds of ratios and pleas for strict 
limitation of its meaning. Like any group of lexicogra- 
phers, the subcommittee has tried neither to walk too far 
ahead of usage nor to lag too far behind. It is apparent, 
however, that a more positive approach is needed. 
There have been several proposals for taking care of 
this situation. M. W. Baldwin, Jr. and R. E. Graham, in 
an unpublished memorandum dated March 24, 1947, 
suggested that the name “decilu” (derived from decilo- 
grathmic unit), with abbreviation dl, be employed to 
express the ratio between two values of any quantity, 
the number of decilus being equal to 10 times the com- 
mon logarithm of the ratio. While many persons sub- 
scribed in principle to this proposal, the decilu has not 
caught on. Subsequently J. W. Horton, in an excellent 
treatment of this problem suggested that the name 
“logit” (an abbreviation of logarithmic unit) be used in 
the same way, i.e. the number of logits being equal to 
10 times the logarithm of the ratio of two magnitudes of 
any quantity. The suggested abbreviation for the logit 
was “Igt.” Despite the fact that this proposal received 
strong encouragement from the Standards Committee of 
the I.R.E., the name logit has not taken hold. Others 
have proposed that the “brigg” be defined as equal to 
the common logarithm of any magnitude ratio, with the 
“decibrigg” equal to 10 times the common logarithm of 
the ratio. Still others have suggested that the name 
“decade” be employed in the same sense as proposed 
above for brigg. These latter proposals have met the 
same fate as their predecessors. The best available 
terminology today would seem to be the rather cumber- 
some and imprecise expression that two values which bear 
a ratio of about 10 to 1 differ by an order of magnitude. 
There seems to be rather widespread agreement that 
a new logarithmic unit for expressing the ratio of two 
values of any kind of quantity is needed. Many people 
would agree, too, that the preferred magnitude for the 
new unit is 10°-! or, in other words, that the number of 
units should equal 10 times the common logarithm of 
the ratio. Why then have none of the proposals to date 
gained currency? The answer seems to be that defining 
bodies are reluctant to accept a new term until it has 
achieved fairly general usage, while in the case of any- 
thing as fundamental as a unit of measurement, most 
people want the stamp of authority before accepting 
something new. The solution obviously lies in the same 
procedure which has been used in the past for many 
kinds of units, namely, official authentication by na- 
tional and preferably international authorities as a war- 


- rant for general use. While defining bodies are usually 
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best advised to let usage develop naturally, new units 

are an important exception since early standardizing 

action can prevent imprecise and inconsistent usage. 
It seems clear, then, that a new unit is needed to 
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represent the ratio of two values of any quantity. The 
size of the unit should be 10°-1, which is to say that the 
number of units should equal 10 times the common loga- 
rithm of the ratio. There remains the choice of a name 
for the new unit. The name “logit,” which has enjoyed 
the greatest publicity, is open to certain objections. If 
pronounced with a hard g it is not euphonious, nor is it 
etymologically correct since logit logically should have a 
soft g. In addition, the proposed three-letter abbrevia- 
tion “Igt” is less desirable for both spoken and written 
use than a two-letter abbreviation analogous to db. The 
name “decilu” is good but the abbreviation “dl” is 
awkward since the same symbol on the typwriter de- 
notes a lower case letter | and the numeral 1. Brigg would 
bring with it the decibrigg, together with its abbrevia- 
tion “db” which is untenable because the same abbrevi- 
ation stands for decibel. Decade, which is inconsistent 
with the present accepted usage of this word, would 
be a most unfortunate choice. 

In view of the nonacceptance of names thus far pro- 
posed, it seems better to try a new one. The name 
should be reasonably euphonious and should lend itself 
to a simple and unmistakable two-letter abbreviation. 
Search of various possibilities leads to the suggestion of 
the name “decilog” (meaning one tenth of the loga- 
rithm) and the abbreviation “dg.”! The similarity of the 
name and the abbreviation to decibel and db, respec- 
tively, should facilitate understanding and acceptance. 
The use of the initial and final letters to form an ab- 
breviation is analogous to stock ticker practice where 
SY stands for Sperry Corporation, AN for Air Reduc- 
tion Co., and so forth. The theoretical and practical con- 
siderations which have been so ably advanced by J. W. 
Horton in his article on the use of relative magnitudes 
are all applicable to this proposed new unit. 

Closely related to the matter of a logarithmic unit for 
expressing magnitude ratios is the need for a logarithmic 
way of expressing the ratio of a magnitude to some 
standard magnitude of the same quantity, or in other 
words, expressing the level of a quantity with respect 
to some reference level. In the case of the decibel certain 
abbreviations have become current for this purpose. 
Thus dbm expresses the power level in decibels referred 
to 1 milliwatt and dbw the power level in decibels re- 
ferred to 1 watt. While such abbreviations are con- 
venient, the application of the new term “decilog” to 
many different quantities will make it difficult to find 
suitable abbreviations. It is proposed, therefore, that 
levels be designated by writing the reference level in 
parentheses immediately after the abbreviation dg. 
Thus 40 db (1 gauss) would signify a level of magnetic 
flux density 40 dg above 1 gauss, 1.e. 10,000 gausses. 

In summary, it is proposed to request the approval of 
standardizing bodies for the definition of a new term as 
follows: 

1 Since writing the above, the writer has learned that A. G. Fox 


in 1951 proposed the same unit, to be named “decilog” with abbrevia- 
tion “dl.” This coincidence suggests that the name is a happy choice. 
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DECILOG (dg) 

The decilog, abbreviated dg, is a dimensionless unit 
for expressing the ratio of two values of any quantity, 
the number of decilogs being 10 times the logarithm 
to base 10 of the ratio. 

The magnitude or level of a quantity may be stated 
in dg, expressing the ratio to a reference value. Thus 
N dg (1 gauss) signifies a level of magnetic flux den- 
sity NV dg above 1 gauss. 


PGAP 


As of December 31, 1953 there were 1,219 paid mem- 
bers of the Group. After a canvass, all unpaid members 
have been removed from the enrollment list. 

Welcome to our newest chapter: Albuquerque—Los 
Alamos! R. B. Jacques is Chairman. 


LocaAL CHAPTERS 


The Chicago Section of the PGAP, under the general 
Chairmanship of G. P. Kearse and Papers Chairman 
L. Krahe, have had an active year of meetings. The fol- 
lowing meeting and talks have been held in the past 
months: 

Oct. 16, 1953—“Graphical Integration of Patterns for 
Antenna Gain Measurements,” by Dr. R. R. H. 
Yang, Andrew Corporation. 

Dr. Yang discussed the basic theory involved in 
determining antenna gain by pattern integration 
methods. The mathematical formulas were derived 
and used in the construction of a special type of 
graph paper. Dr. Yang described the method of 
using this paper to eliminate many of the math- 
ematical computations involved. A question and 
answer period followed. 

Nov. 20, 1953—“UHEF TV Receiving, Antenna Meas- 
uring Techniques,” by Dr. R. M. Soria, American 
Phenolic Corporation. 

Dr. Soria presented a slide film showing propaga- 
tion characteristics of vhf and uhf TV signals and 
comparing the performance characteristics of 
standard vhf and uhf TV receiving antennas. The 
film was followed by a brief description of vhf and 
uhf test equipment and measuring techniques. A 
question and answer period followed. 

Dec. 18, 1953—(Joint Meeting with Microwave The- 
ory and Technique Group)—‘“Waveguide Coupled 
Circuits,” by Mr. Tore N. Anderson, Airtron, Inc. 

Mr. Anderson covered the theoretical and ex- 
perimental aspects of all basic types of directional 
couplers including the Bethe Hole Couplers, Mo- 
reno Couplers, Crossed Guide Couplers, and Multi- 
Hole Couplers. Mr. Anderson explained the limi- 
tations, uses, relative advantages, and gave demon- 
strations of each type. 

March 19, 1954—“UHF TV Propagation Problems,” 
by Mr. H.N. Frihart, Motorola. 
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MEETINGS 


The Los Angeles Chapter of the PGAP is organizing 
three sessions of technical papers on antennas and prop- 
agation at the August 1954 Western Electronics Con- 
vention in Los Angeles. 

Our Professional Group is sponsoring a Spring Tech- 
nical Meeting jointly with the USA National Committee 
of the International Scientific Radio Union (URSI). 
This meeting will be held at the National Bureau of 
Standards, Washington, D. C. on May 3, 4, 5, and 6, 
1954. 

A meeting of the USA National Committee and meet- 
ings of participating commissions will be held on Mon- 
day, May 3. The technical sessions will be held on Tues- 
day, Wednesday, and Thursday, May 4, 5, and 6: 

The participating Commissions and their chairmen 
are: 

Commission 1: Radio Measurement Methods and 

Standards 
Dr. Rufus G. Fellers, 

Bldg. 3, PH, Radio 1, 
Naval Research Laboratory, 
Washington 20, D. C. 
Commission 2: Tropospheric Radio Propagation 
Chairman: Dr. Henry G. Booker, 
School of Electrical Engineering, 
Cornell University 
Ithaca N.Y: 
Commission 3: Ionospheric Radio Propagation 
Chairman: Dr. L. V. Berkner, - 
Associated Universities, Inc., Room 
6920, 
350 Fifth Ave., 
New York 1, N. Y. 
Commission 4: Torrestrial Radio Noise 


Chairman: 


Chairman: Mr. Frederic H. Dickson, 
Communications Liaison Branch 
OCSigO, 
Room 2D276, National: Defense 
Building, 


Washington 25, D.C. 

Radio Astronomy 

Dr. John P. Hagan, 

Naval Research Laboratory, 
Washington 20, D. C. 

Radio Waves and Circuits, Including 
General Theory and Antennas 

Profi. C.. Jordan, 

Department of Electrical Engineer- 

ing, 

University of Illinois, 

Urbana, Ill. 

Electronics 

Prof. W. G. Shapherd, 

Institute of Technology, 

l _ University of Minnesota, 

Minneapolis 14, Minn. 


Commission 5: 
Chairman: 


Commission 6: 


Chairman: 


Commission 7: 
Chairman: 


ey 


1954 


The success of these meetings depends upon the 
whole-hearted co-operation of the individuals partici- 
pating in them and a continuation of the enthusiastic 
response we have had in previous years. 

Advance registration forms and details of the meeting 
will have been sent to you on or about April 1, 1954, to- 
gether with the preliminary program of titles. 

Your co-operation in making these URSI-IRE meet- 
ings a success in the past and your continued effort an- 
ticipated in similarly producing a worthwhile meeting 
this Spring are greatly appreciated. It is hoped that you 
and other members of your organization will be able to 
attend and participate in the meetings. 


CORRESPONDENCE 


The editors of the PGAP Transactions have re- 
cently received a letter from W. E. Gordon of Cornell 
University surveying the IRE and the Professional 
Group System in general and the PGAP in particular. 
The letter is reproduced in its entirety for the informa- 
tion of the membership. 

“A professional society should provide its members 
with (a) a meeting place for the presentation and dis- 
cussion of original contributions to the science, (bd) a 
journal for the publication of these contributions, and 
(c) a publication devoted to reviews and news. 

“The Institute of Radio Engineers because of the di- 
versified interests of its 37,000 members is not a profes- 
sional society but should be considered as a collection of 
professional societies. The professional group movement 
is evidence of this. The professional groups must assume 
the status of professional societies with strong assistance 
from the Institute. The assistance should be substantial 
judging from the current financial statement of the In- 
stitute. 

“If the PGAP is to become a professional society it 
must fulfill the three requirements set forth in the first 
paragraph. PGAP is providing the meeting place 
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through joint sponsorship of URSI meetings, although 
membership in PGAP is not required of persons attend- 
ing the joint meetings. The TRANSACTIONS OF THE IRE 
PGAP does not constitute the journal and, it seems to 
me cannot, for two reasons. First, the scientific scope is 
limited. No member of PGAP limits his radio interests 
to the narrow field of antennas and propagation, and 
hence the TRANSACTIONS do not cover the interests of 
the membership. The scope should be broadened to in- 
clude what might be called radio physics embracing ra- 
dio waves in the atmosphere, both troposphere and up- 
per atmosphere, the physical problems associated with 
radio waves in the atmosphere, and radio waves from 
extra-terrestrial sources. I realize that to broaden the 
scope of the TRANSACTIONS requires broadening the 
scope of PGAP, but feel that this is required not only for 
the reason just cited but, secondly to enlarge the mem- 
bership of the Group so that a first class journal will be 
economically feasible. 

“The PROCEEDINGS could serve a useful purpose to the 
Group and all other members by expanding its review 
and tutorial paper policy, and disseminating news of 
interest, thus fulfilling requirement (c). 

“In becoming a journal, the TRANSACTIONS OF THE IN- 
STITUTE OF RADIO ENGINEERS PROFESSIONAL GROUP 
ON ANTENNAS AND PROPAGATION must change its name 
to something less awkward and something that will be 
referenced more usefully in libraries. Radio Physics or 
Radio Waves seems more appropriate. 

“Since the TRANSACTIONS is just starting I feel strongly 
that now is the time to make the change in a direction 
which more satisfactorily serves the interests of the 
members and which promotes survival both on a scien- 
tific and economic basis. To this end the editor with the 
approval of the Professional Group should formulate a 
policy regarding the scope and format and specifically 
invite papers representative of the new scope of the 
journal and the Professional Group should secure the 
support of the Institute.” 


OC) 
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contributions 


- An Experimental Investigation of the Single- 
Wire Transmission Line* 


T. E. ROBERTS, Jr.f, ASSOCIATE, IRE 


Summary—The results of some measurements made on a di- 
electric coated wire are presented and compared with theoretical 
results. These measurements indicate that the single-wire line can 
be considered as a simple transmission line provided account is taken 
of the ‘‘end-effect.” 


INTRODUCTION 


INCE THE INTRODUCTION of the dielectric 
S coated wire as a practical means of transmitting 

high-frequency radio waves,! there has been con- 
siderable speculation as to whether transmission-line 
concepts were really applicable to such a structure. It 
has recently been shown?? that the structure can be 
treated as a transmission line whenever its length is 
sufficiently great. 

The purposes of the present research have been to 
perform some rather fundamental experiments that 
clarify the theoretical results, and also to develop meas- 
uring techniques applicable to the single-wire line. 

The physical quantities of most interest are: 

1. attenuation and phase velocity of the principal 

mode, 

2. equivalent circuits for the transition from a co- 

axial line to the single-wire line, 


* Manuscript received by PGAP, March 22, 1953. Part of a thesis 
submitted in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy. This research was made possible by support 
extended Cruft Laboratory, Harvard University, under contract 
N5ori-76, T.O. 1. 

} Formerly, Harvard University; now, American Machine and 
Foundry Co., Greenwich, Conn. 

1G. Goubau, “Surface waves and their application to transmis- 
sion lines,” Jour. Appl. Phys. vol. 21, p. 1119; 1950. 

2 G. Goubau, “On the excitation of surface waves,” Proc. I.R.E., 
vol. 40, p. 865; 1952. 

3 T. E. Roberts, Jr., “Theory of the single-wire transmission line,” 
Jour. Appl. Phys., vol. 24, p. 57; 1953. 


3. current and charge distributions, 

4. radiation-field patterns. 

The first three quantities can all be measured with 
the same equipment. Measurement of the field pattern, 
however, requires entirely different equipment, and in- 
volves a considerable project. For this reason, and be- 
cause our primary interest is in the transmission quali- 
ties of the single-wire line, it was decided not to make 
measurements of the field patterns. 
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WIRE WITH A SURFACE MODIFICATION 


RF 


GENERATOR; sain ag 


DETECTOR 


Fig. 1—Block diagram of the system. 


The required equipment is quite simple. It consists of 
a stable radio-frequency generator,* detecting equip- 
ment, and a straight wire supported at each end, with 
a traveling probe or loop to measure the magnitudes of 
the electric and magnetic fields. A block diagram of the 
system is shown in Fig. 1, and a photograph of the ex- 
perimental arrangement is shown in Fig. 2. 

A wavelength of 3.2 cm was chosen for the radio- 
frequency signal. This appeared to effect a reasonable 


4 T. E. Roberts, Jr., “A stable amplitude-modulated microwave 
generator,” Proc. I.R.E., vol. 40, p. 1125; 1952. 
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Fig. 2—The experimental arrangement, 


compromise between accuracy and over-all size, making 
it possible to obtain a probe-travel over the entire length 
of line (about 80A). An over-all length of 80d is sufficient 
for making measurements on wires that are of practical 
interest, although all properties of unmodified wires 
(excepting wires of small radius) cannot be measured. 

Accurate methods of measuring charge and current 
distributions on cylindrical structures, developed by 
Morita,® make use of small probes and shielded loops 
extending from a slotted tube. This method cannot con- 
veniently be applied in the present case because of the 
small dimensions involved. Also, the use of long slotted 
tubes would greatly reduce the possible choice of wires 
that could be investigated. It was therefore decided to 
use probes and loops that were supported from the out- 
side by the smallest coaxial line available (1/32 inch 
o.d.) and to accept the consequent errors. 

A photograph of the probe and shielded loop is shown 
in Fig. 3. It will be noticed that neither of these units is 
really small compared to the wavelength, but it can be 
shown that the error introduced by their size is of con- 
sequence only near the ends of the wire. A rather serious 
error of another type, commonly referred to as “probe- 
loading,” is introduced by the supporting coaxial line 
and the loop (or probe) considered as a unit. This unit 
affects the charge and current on the line—that is, it 

5 T. Morita, “Measurement of Current and Charge Distributions 


on Cylindrical Antennas,” Cruft Laboratory Technical Report No. 
66, Harvard University; February, 1949. 


acts as a load (predominantly reactive in this case). 
This second error cannot be completely eliminated in the 
measurement of current and charge, although it can be 
reduced somewhat by using loose coupling. If the loop 
is spaced any great distance from the wire there will be 
a considerable error introduced near the ends. In the 
case of a long wire there can be a rather large spacing 
between the measuring loop (or probe) and the wire, 
without introducing an appreciable error in the dis- 
tributions except near the ends. In this way the “probe- 
loading” effect is made negligible for measurements 
made away from the ends. In making measurements 
near the ends the coupling is increased, and there may be 
some error from “probe-loading.” Measurements can 
be made by making use of this “probe-loading” effect 
in a scheme developed by Goubau.* The method is sim- 
ilar to the method of measuring field strengths in 
resonant cavities. It is not used in this work, however. 

The purpose of the circular plate on the charge probe 
(Fig. 3, p. 48) is to isolate the probe from the supporting 
structure. Without such a plate the response of the 
probe was inconsistent and a number of errors were 
introduced. The small cap on the tip of the probe in- 
creases the coupling to the wire and appears to give a 
slightly better measure of the radial electric field close 
to the wire. The loop pick-up is shielded and symmetri- 
cal; a plate below the loop would serve no use. 


6 G. Goubau, “Zur ausmessung elektromagnetischer felder nuttels 
testkérper,” Hochfreq. und Elektroak., vol. 62, p. 73; 1943. 
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Fig. 3—The current loop and charge probe. 


The probe carriage is mounted on a steel track by 
means of four flanged (railroad-type) wheels. The posi- 
tion of the probe is indicated in centimeters by a steel 
tape. For accurate distance measurements over short 
intervals a vernier is mounted on a small magnet that 
may be placed anywhere along the steel track without 
danger of slipping. The carriage is shown with the other 
equipment in Fig. 2. 

In some of the measurements coaxial horns were used 
to make the transitions from the single-wire to the co- 
axial lines. Two identical horns were constructed from 
sections of 30-degree cones of 0.003-inch brass stock. 
The cones tapered out from a diameter of 0.556 inch to 
a diameter of 3.1 inches. At the point where the cone 
joins the coaxial line of the same diameter there is a 
sharp junction which probably gives considerable re- 
flection. 


SUMMARY OF THEORETICAL RESULTS 


The specific structure of interest here (see Fig. 4) 
consists of a finite length of wire (with a modified sur- 
face) connected to a flanged coaxial line at each end. A 
number of theoretical results for this structure were re- 
ported earlier, and a brief summary of them is given 
below. 


April 


The total axial current in the wire J(z) may be repre- 
sented as a sum of two components, Jo(z)+J’(z). Lo() 
is a “transmission-line current” with a propagation con- 
stant fo> (6 is the free-space propagation constant). 
The auxiliary component of current J’(z) is a compli- 
cated function of z. In general, however, some simpli- 
fication occurs if the length of the line L is sufficiently 
large. In this case the current J’(z) is of significance only 
near the ends of the wire. This auxiliary current is quite 
large at each flange and decreases to small values at 
points some distance from the flange. Theauxiliary cur- 
rent is a result of the voltage V in the gap and is pro- 
portional to this voltage. For the case where the gap is 
small, numerical values of the component of current 
I(z) were computed for two specific cases* and these re- 
sults used to compare with the experimental values in a 
later section. 

It was also shown that the coupling between the two 
coaxial lines was by means of the principal mode,’ that 
is, the equivalent circuit is that shown in Fig. 4. The 
parameter G; is the radiation conductance; efficiency of 
coupling is directly related to the ratio Gi/Go. 

The apparent admittance terminating the single-wire 
line is indicated by the standing-wave ratio, which must 
be measured at some distance from the ends of the line 
in order to be a measure of the propagating mode only. 
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Fig. 4—The single wire connected between flanged coaxial lines. 


Since a coaxial-horn feed introduces a smaller auxiliary 
current than the sharp transition of the flange, this 
type of feed is often used at the transmitter end to ex- 
cite the line. If the voltage at the receiving end is zero 
(i.e., a short circuit), no auxiliary current is introduced 
at that end. Then the attenuation and phase velocity 
of the wire can be determined by measuring the parame- 
ters A, +jF,=jfostp.+j®,, at a distance s from the 
short-circuited end. {os represents the attenuation and 
phase shift of the line; p, and ®, are the terminal 
functions. 
7 This is true only for large L. 


8 Ronald King, “Transmission-line theory and its application,” 
Jour. Appl. Phys., vol. 14, p.577;1943,, 
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PHYSICAL PROPERTIES OF THE EXPERIMENTAL LINES 
Most of the measurements presented here were made 
on two dielectric-coated wires, No. 1 and No. 2, having 
the characteristics given in Table I. A few measure- 
ments were made on wire No. 3, also described here. 
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may be in error in this case. Before discussing the meth- 
ods used in arriving at the experimental values of Table 
II, it is convenient to discuss some of the other measure- 
ments. This discussion leads to the procedures used to 
obtain the above data. 


TABLE I 
PHYSICAL PROPERTIES OF DIELECTRIC CoATED WIRES 


Wire No. 2 Wire No. 3 
0.0572 cm a’ 0.0318 cm 
0.0610 cm a 0.0279 cm 


Wire No. 1 
a’ 0.159 cm a’ 
a 0.203cm _ a 
ex/e 225 \ i 
ian} 0.0010 ;measured ns 
o 1.5107 mhos (Brass) o 
Dielectric Material Vinyl Plastic “Spaghetti” 


—, 


2.50 
a, (M.I.T. Tables) 
5 


.65 107 mhos (hard-drawn copper) 


Dielectric Material—Formvar (a trade name, General Electric Company) 


The properties of the dielectric coating of wire No. 1 
were measured at the operating wavelength of 3.2 cm 
by placing a short piece of hollow “spaghetti” coating 
in a rectangular cavity and observing the shift in the 
resonant frequency and the change in the Q of the cav- 
ity. From these measurements the values of the relative 
dielectric constant, ¢€;/€), and the loss tangent, tan 6, 
are easily calculated. 

Since the dielectric coating on wire No. 2 is very thin 
and not readily removed, it would be difficult to meas- 
ure accurately its dielectric constant. Previous meas- 
urements have been made on a block of this material, 
which has the trade name of Formvar, by the M.I.T. 
Laboratory for Insulation Research, and the data given 
above were taken from their report.® It is not certain 
that the dielectric coating used here is identical with the 
sample tested at M.I.T., but no other data on this ma- 
terial are available. If there is any slight variation in 
the materials, one usually expects the loss tangent to 
be quite different and the dielectric constant to be rela- 
tively unchanged. 


CURRENT AND CHARGE DISTRIBUTIONS 


It has been shown elsewhere’ that a finite wire of 
sufficient length can be terminated so that it behaves 
as a flat line except near the ends. Under terminated 
conditions the current and charge distribution near the 
driving end of this finite line should be the same as for 
the infinitely long wire. 

An experiment was arranged with the transmitter 
connected as shown in the insert of Fig. 5. A horn was 
connected at the receiving end, and the load tuned until 
a minimum standing wave was indicated by the travel- 
ing probe near the center of the line. The probe was 
then moved to the receiving end and the line appeared 
almost flat (swr <1.05) except in a small region near 
the receiving horn. Near the generator the current rises 
sharply as expected from the theory. The result of the 
measurement on wire No. 1 is shown in Fig. 5 together 
with a theoretical curve. (The theoretical curve was ob- 
tained by adding the two components of current J (z) 
and I’(z). Curves of I’(z) for wires No. 1 and No. 2 have 
been given elsewhere. ?°) 


TABLE II 
TRANSMISSION LINE PROPERTIES OF DIELECTRIC COATED WIRES 


wire rola co af nike Gi/Ge 
No. 1 0.092 0.40 0.0079 0.0062 0. rg 
No. 2 0.010 0.12 0.0039 0.0040 “Be 


a die a cu a Total Not 

db/m db/m db/m cm 

0.43 0.17 0.60 3,120 Theoretical 
0.53 Shesilat Experimental 

0.06 0.125 0.185 3.185 Theoretical 
0.24 ails) Experimental 


ooo 


¢ Guide wavelength. 


Table II contains a summary of the most basic data 
related to these two wires. The experimental values for 
both wires are seen to be in rather good agreement with 
the theory. One should expect somewhat closer agree- 
ment between the measured and theoretical values of 
the attenuation of wire No. 2, but, as mentioned above, 
the value of the loss tangent of the dielectric coating 


® “Tables of Dielectric Material,” Report V, M.I.T. Laboratory 


for Insulation Research; February, 1944. 


The theoretical current-distribution curve of wire No. 
1 is interesting in that there is a slightly dip below the 
final constant value. This dip can be interpreted as - 
interference between the two components of current J . 
and J’, since J) has a propagation constant of {> and I’ 
has an “average propatation constant” of 8, the free- 
space value. 


10 Roberts, “Theory of the single-wire transmission line,” <bid., 
Figs. 5 and 6. 
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Fig. 5—The current on a terminated single wire, wire No. 1. 


The accuracy of the measured curve of Fig. 5 is 
rather poor because the result depends directly on the 
probe-to-line spacing. Moreover, some error is intro- 
duced by the small residual standing wave and by probe 
loading. Considerable improvement was made by sup- 
porting the line with polyfoam so that the unsupported 
span was only about 50 cm. Even so, it cannot be con- 
cluded from the measurements that the dip exists. 

The lack of accuracy can be reduced by measuring 
the distribution in a somewhat indirect manner. Let 
the positions of the generator and load be interchanged. 
The current near the receiving end (flanges in this case) 
is closely related to the current near the generator on an 
infinitely long wire, as is seen from the following. 

If a load, Yz, is connected to terminals 11 such that 
Y,+Yi1=Gp, there will be no reflected transmission- 
line current, in other words, the wire is terminated. 
There will be a standing wave near the flange caused 
by a current of type I’(z). The total current will be 
I(z) =Iye—#% + I'(z). Except for the reversal of the prop- 
agation direction of the current [(z), this expression is 
identical to the current on an infinite line. The two 
components of current interfere and produce a standing 
wave, as shown in Figs. 6 and 7. The agreement be- 
tween experiment and theory is seen to be quite good. 
The maximum discrepancy occurs for small values of z. 
The theoretical curves are expected to be more accurate 
than the experimental values because of possible probe- 
loading effects. It seems, however, that the relative mag- 
nitude of I (z) can be determined to a reasonable accu- 
racy from such an experiment. Any variation in probe 
coupling can be easily corrected because the current 
distribution must be symmetrical about a horizontal 
line:! 


The curve of Fig. 6 does not show a slow interference 
beat as occurs in the theoretical curve of Fig. 5 because 
there is no reflected transmission-line current I o(z) 
= I e%o?. Such a beat is detected if the load is tuned to 
some value other than Y;=Gp— ¥;. The distribution on 
wire No. 1 for Yr = —7B, is shown in Fig. 8, p. 52, where 
interference between the reflected components of 
current [ye%o? and I’(z) is quite evident. The measured 
distributions on wire No. 1 fora larger gap, b=0.22, 
with Y,;=jB are given in Fig. 9, p. 53. Note that the 
final standing-wave ratio is much greater for the large 
than for the small gap. This is to be expected since the 
radiation conductance, G, is less for the large gap; that 
is, there is a stronger coupling between the coaxial line 
and the single-wire line. 

The above measurements are in quantitative agree- 
ment with the theory and are a good illustration of the 
advantage gained by separating the current into two 
components. 


MEASUREMENT OF ATTENUATION AND 
PHASE VELOCITY 


The experimental measurements given above verify 
the theoretical result that the single-wire line, if it is 
sufficiently long, may be treated.as a transmission line. 
An idea of its required length has also been obtained 
in the two specific cases considered. It would appear, 
therefore, that the transmission-line techniques of 
measurement may be applied provided care is taken 
in application that ranges of validity are not exceeded. 


1 In the measurements presented here a correction curve was de- 
termined by short-circuiting the gap and measuring the magnitudes 
of the maximum of the resulting high standing wave. From these data 
the probe was calibrated. SE ri 
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Fig. 6—Current and charge distribution, wire No. 1, Yr =Go— Y; (line terminated); b=a. 


8 9 10 aes ees 


X CHARGE 
e CURRENT 


} EXPERIMENTAL 
CURRENT -THEORETICAL 


I,(0) = 3.94 MO. 


VOLT 


Fig. 7—Current and charge distribution, wire No. 2, Y,=Go—Y; (line terminated) b =a. 


One of the most accurate methods of measuring small 
attenuation is to measure the standing-wave ratio on a 
short-circuited line. In the present work, the wire was 
terminated in a large rectangular flange (11.2 X14.4n) 
which was expected to represent a good short circuit, 


and a coaxial horn was used to couple to the transmit- 
ter. The standing-wave ratio was determined at points 
about two meters from the short circuit by measuring 
the “width” of the minimum. The relation between at- 
tenuation and the above “width measurements” is 
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Fig. 8—Current and charge distribution, wire No. 1, Y= —jB; b=a. 


given by 
VJ p? — 1 sinh as = sin £6, 


where # is the ratio of the current (charge) at a distance 
6 from the -‘minimum to the current (charge) at the 
minimum. Since a horn was used to couple into the line 
and the measuring probe was located near the center of 
the line, the current J’(z) may be considered negligible, 
and the results should be an accurate measure of the 
attenuation. The guide wavelength (equivalent to phase 
velocity) was determined from the distance between 
minima. These results are given in Table II. 

Another method of measuring the attenuation yielded 
results that are substantially in agreement with the 
measurements just discussed. Here, two horns are con- 
nected to the line and the received signal is measured 
as a function of the horn separation, L. The final slopes 
of the curves of received signal vs separation should be 
equal to the attenuation. These curves are plotted in 
Figs. 9 and 10 and the attenuation measured by the 
two methods is indicated. Some deviation from a 
straight line for small spacing is of course to be expected 
since the fields are quite complicated near each horn. 

It is interesting to compare the data of wires No. 1 
and No. 2 with data taken for similar wires with no 
dielectric coating. It is seen in Figs. 9 and 10 that in the 
case of no dielectric coating the horns are not sufficiently 
far apart for the attenuation of the received signal to 
approach a simple linear function of the distance. The 
measured curve is expected to approach an asymptote 
of the slope indicated by the theoretical attenuation 
curves. The length of wire is not sufficiently long to 
show this. Thus it can be stated that.the uncoated wires 
of the length used here cannot be classified as transmis- 
sion lines but the coated wires can be so classified. This 
comparison of the coated and uncoated wires clearly 
displays the guiding properties of a thin dielectric coat- 
ing. (The dielectric coating was only 0.015 inch thick 
for wire No. 2.) 


THE RADIATION CONDUCTANCE Gj 


The transition from a flanged coaxial line to a single- 
wire line may be represented by a four-terminal net- 
work. If the gap is sufficiently small, so that the radial 
electric and magnetic field distributions in the gap are 
the same whether the circuit is driven from the coaxial 
line or from the single-wire line, then the general four- 
terminal network reduces to a simple shunt circuit. 
This shunt circuit is located at the flange and may be 
represented by an admittance Yi=Gi+jBi plus an 
auxiliary transformer. 

Suppose the network is represented by a simple shunt 
circuit. The most important parameter is Gi, since the 
susceptance B and the transformation ratio m can be 
changed by tuning adjustments made inside the coaxial 
line. These tuning adjustments (by use of lossless ele- 
ments) do not change the value of G;. The relative con- 
ductance G;/Go (where Gp is the characteristic conduct- 
ance of the single-wire line) is a measure of the efficiency 
of coupling between the two lines.?? 

The radiation conductances (G;) of wires No. 1 and 
No. 2 were determined by the Weissfloch¥- method 
by placing a movable short circuit inside the coaxial 
line and measuring the input admittance to the single- 
wire line at a point about two meters from the junction. 
(The radio-frequency excitation, \ = 3.2 cm, was coupled 
into the single-wire line by means of a horn.) The meas- 
ured input admittance was transformed to the reference 
plane at the junction of the two lines by subtracting 
out the attenuation of the length of line used. The re- 

” For example, the efficiency of coupling (on a power basis) from 
a coaxial line to a terminated single-wire line is 1/(1+G,/Go), and the 
efficiency of coupling from a single-wire line to a coaxial line which 
terminates the single-wire line is (1—G,/Go). 

*% A. Weissfloch, “Ein transformation tiber verlustose vierpole 
und seine anwendung,” Hochfreq. Electroak., vol. 60; 1942. 

4 N. Marcuvitz, “On the representation and measurement of 
waveguide discontinuities,” Proc. I.R.E., vol. 36, p. 728; June, 1948. 

% A. A. Oliner and H. Kurss, “The Precision Measurements of 


the Equivalent Circuit Parameters of Dissipative Microwave Struc- 
reve Paper No. 78, I.R.E. National Convention, New York, N. Y.: 


a 
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Fig. 9—Received signal as a function of separation, L, wire no. 1. 
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Fig. 10.—Received signal as a function of separation, L, wire no. 2. 


sulting input-admittance measurements were then 
plotted on the Smith chart and the parameters G, and G, 
were determined from the best circle. The value of G 
was zero (within the experimental accuracy) for all di- 
ameters of the coaxial line used. G,, therefore, repre- 
sents G, directly. An example of the data for wire no. 1 
is given in Fig. 11, p. 54. 


In this way values of G; were determined for a num- 
ber of gap radii, b; these results are in Fig. 12, p. 54. 
A curve is drawn through the measured points, then 
extrapolated to the value at zero gap, which value it 
approaches at zero slope. 

Measurements made for very small gaps are not ac- 
curate because of the losses in the coaxial line. The ef- 
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Fig. 11 (above) —Experimen- ee rie 28, 
tal admittance plot wire No. 1; ne ee 
b/a =2.45, 
fect of these losses is easily determined by changing 
the length of the coaxial line by several half-wavelengths ra : Ss Septet ee 

. . 0. = 
and repeating the measurement. If the losses in the co- a 


axial line are appreciable, the measured values will be ; ieieud eae 
different in the two cases.!® These losses were quite Ge bes sans 
noticeable for gaps of about half the smallest size indi- BS 
cated in Fig. 12. However, no such “loading effect” was Be 
detectable for the values given in these figures. 

The results of the above extrapolation are tabulated bs 


with the theoretical values in Table II. The agreement 
is very close. It is, in fact, somewhat better than is ex- 
pected from an inspection of the experimental data of 
Fig. 11. 


_ * Montgomery, Dicke, and Purcell, “Principles of Microwave 
Circuits,” M.I.T. Radiation Laboratory Series, vol. 8, McGraw-Hill 
Book Co., New York, N. Y., p. 231; 1948. 
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Fig. 13—Data for determining the insertion loss of the coaxial horn; wire no. 1. 


THE MEASURED EFFICIENCY OF THE COAXIAL HORNS 


As a point of practical interest, the efficiency of the 
coaxial horn was determined experimentally. The horn 
was connected to wire no. 1 at the receiving end in place 
of the flange, and the measurement was carried out as 
in the previous section. The data are given in Fig. 13. 
The results of the measurement are: insertion loss, 0.17 
db—or efficiency, 96 per cent. 

Since, for this specific case, the center of the measured 
admittance circle lies at the center of the Smith chart, 
the equivalent circuit can be represented by an atten- 
 yation of 0.17 db, plus a properly located transformer. 


After carrying out the complete procedure the trans- 
formation ratio was found to be »=1.05. Thus the 30- 
degree, 6-inch-diameter coaxial horn is a reasonably 
efficient device with very low insertion and reflection 
losses. 

The measurement was also made in the reverse sense, 
i.e., by connecting the short circuit, a large conducting 
plate, to the single-wire line and measuring the admit- 
tance seen by the coaxial line. Of course, the measured 
admittance circle was not centered at the center of the 
Smith chart, but the data were in agreement with the 
equivalent circuit found above. 

In general, more accurate data are obtained by mak- 
ing the measurements in the reverse sense noted above, 
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although there is additional graphical work involved in ACKNOWLEDGMENTS 
subtracting out the attenuation of the single-wire line 
in this case.'” 


The author wishes to acknowledge the encourage- 
ment and guidance of Professor R. W. P. King and the 
valuable suggestions of Dr. T. Morita and others of the 
Cruft Laboratory. 


17 The single-wire line must be rather long (usually causing ap- 
preciable attenuation) in either case, in order to avoid the auxiliary 
current effects. 


Sweep Frequency Backscatter—Some 
Observations and Deductions” 
RICHARD SILBERSTEIN{ 


Summary—Sweep-frequency backscatter records have proved to 
be of great value in identifying the sources of backscatter seen on a 
fixed frequency by demonstrating the development of the echo as 
frequency and range increase. The most commonly observed scatter 
is ground scatter propagated via the F2 layer, but it is also evident 
that the other layers propagate ground scatter and that scatter from 
the distant E region may at times be important. 

In one group of observations over an 1,150-km path on three 
undisturbed days, the values of F2-layer maximum usable frequency 
scaled from midpoint vertical-incidence ionospheric records and 
those determined by backscatter delay assuming ground scatter 
agreed almost within experimental error. In another three-day group 
characterized by a low-latitude ionospheric disturbance with low 
geomagnetic K indexes but considerable sporadic E activity, values 
of muf determined from scatter were much too high under the 
ground-scatter assumption, errors of about 30% being not un- 
common. 


INTRODUCTION 


SE of sweep-frequency equipment to observe 
characteristics of scatter energy received back at 


the site of a high-frequency radio transmitter rep- 
resents a new approach to the study of the ionosphere. 

Early in 1951 a Model C3 ionosphere recorder! was 
modified for a 3- to 25-mc linear 12-minute sweep and 
connected to two rhombic antennas at the Sterling, 
Virginia laboratory of the National Bureau of Stand- 
ards. These antennas were pointed westward for ob- 
servation of oblique-incidence pulse transmissions over 
the 1,150-km path from Sterling to a site near St. Louis, 
Missouri. An ionosphere recorder was installed at Ba- 
tavia, Ohio, the path midpoint, for obtaining vertical- 
incidence data from which calculations of maximum 
usable frequency and travel time were to be made and 
checked against oblique-incidence results.” 

As early as July, 1951, photographs were made by 
E. E. Ferguson on the Sterling equipment showing back- 
scatter starting locally and running out to about 3,500 
km. Results then were rather poor because of the low 
power, about 30 kw peak, with a 50-microsecond pulse, 
and a receiver bandwidth of 70 kc. 

* Original manuscript received by PGAP, October 5, 1953; revised 
manuscript received January 19, 1954. 


t evar ane i Standards, Washington, D. C. 

. M. Carroll, utomatic ionosphere der,” El } 
vol. 25, pp. 128-131; May, 1952. BP ee et SPREE 
: oe E. pene and P. G. Sulzer, “Sweep-frequency oblique- 
incidence 1onosphere measurements over an 1,150-km path,” Proc. 
I.R.E., vol. 40, p. 1124; September, 1952. i 
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It was later found that, for the specific job of scatter 
studies, a relatively slow rise time and considerable 
pulse elongation could be tolerated, so that a bandwidth 
near the theoretical value for maximum signal-to-noise 
ratio in the presence of white noise (1.2 divided by the 
pulse duration) would be used. Values of 100 usec for 
the pulse duration and 11 kc for the bandwidth were 
finally adopted. In recording the marginal scatter sig- 
nals, base-line shift due to strong interfering signals was 
troublesome. In overcoming this difficulty it was neces- 
sary to apply more video clamping to the oscilloscope 
circuit than was used in the standard Model C3 re- 
corder. Also all short time-constant differentiation was 
eliminated since it contributed to base-line shift and 
also caused gaps following the trace produced by a 
strong echo, which obscured subsequent weak echoes. 
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Fig. 1—Multifrequency sweep for December 15, 
1952, at 1257 EST, showing backscatter. 


Fig. 1 is a photograph of a typical normal-day, sweep- 
frequency back-scatter record taken on December 15, 
1952 at 1257 EST. Perpendicular to the horizontal axis 
are frequency markers at 1-mc intervals running from 3 
mc on the left to 25 mc on the right. The vertical scale 
is height or slant range in hundreds of kilometers, 
running to 4,000 km at the top, ‘the first visible marker 
above the ground pulse at the bottom being for 100 km. 
At the lower left-hand corner is seen the well-known 
vertical-incidence h’ —f trace, which could be obtained 
because some energy does escape from the rhombic an- 
tenna in a vertical direction. Starting at the second- 
order reflection and increasing in range as frequency 
increases is the trace of the scatter echoes presumably 
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coming from the distant ground and propagated via 
the F2 layer.* Starting at 21 mc behind the regular line 
of echoes is an echo at a constant range of 2,500 km. 
Relatively fixed echoes (i.e. echoes not sensitive to fre- 
quency) appear at about the same slant range at the 
higher frequencies on practically all sweeps in which the 
range to the frequency-sensitive scatter is shorter at 
these frequencies. This slant range is correct for the 
Rocky Mountains. The higher frequencies appear to 
work best because of low absorption and favorable low- 
angle antenna characteristics. 

Sweep-frequency “scatter sounding” has several ad- 
vantages over the fixed-frequency type in the study of 
ionospheric phenomena and in the practical evaluation 
of the use of scatter echoes for skip-distance determina- 
tion. By sweeping the frequency it is possible to obtain 
information as to the source or type of scatter by the 
manner in which the echoes “grow” as the frequency is 
increased. Again for the test of any premise, a large 
number of samples are afforded (a) at different fre- 
quencies in any one sweep, and (d) under different di- 
urnal conditions. For instance, in skip-distance de- 
termination over a fixed path there are, in general, only 
two times of day when a fixed frequency skips. With 
sweep-frequency techniques it is likely that there will 
be some frequency which will just skip on the particular 
path at any time of day. 

Present results were limited mainly to daytime hours. 
With the present power and rhombic antenna, the main 
lobes are too high at the relatively low frequencies, 
which have long F2-layer skip ranges at night with low- 
angle requirements. 
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Fig. 2—Multifrequency sweep for December 15, 1952) 
at 1309 EST, with expanded height scale. 


Fig. 2 shows a similar sweep made 12 minutes later 
using a 1,000-km height scale. More detail is visible on 
the manner in which the backscatter develops from the 
second-order F2 echo. The separation in the scatter 
trace which is most prominent at 700 to 800 km and 8.2 


to 9.2 mc seems to be associated with the extraordinary | 


wave.*-4 It is often visible under quiet conditions. The 


3W. Dieminger, “Origin of ionospheric scattering,” Proc. Con- 
ference on Ionospheric Physics, vol. II, 1950, The Pennsylvania State 
College and the Air Force Cambridge Research Laboratories. 
4 Dieminger, “The scattering of radio waves,” Proc. Phys. Soc. B, 


- (London, Eng.) vol. 64, 374B, pp. 142-158; February, 1951. 
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frequency separation of these scatter echoes at near 
vertical incidence is often less than the well-known fre- 
quency separation of ordinary and extraordinary wave 
echoes at vertical incidence. 


A DETAILED StTuDY oF Two DIFFERENT PERIODS 


Analyses were made of records taken during two 
periods which differed radically in ionospheric char- 
acteristics but were both geomagnetically quiet. The 
period January 21-23, 1953 was a quiet period both 
geomagnetically and ionospherically and was marked 
by a very low incidence of sporadic-E reflections. On 
the other hand, the period December 17-19, 1952, 
although quiet geomagnetically, was characterized by 
rapid changes of critical frequency throughout the day 
(characteristic of a low-latitude type disturbance) and 
by a high incidence of sporadic-F reflection. The quiet 
period will be considered first. 


Fig. 3—Sweeps for January 22, 1953. 


Fig. 3 shows a sequence of three sweeps made on 
January 22, 1953. The range scale begins at 200 km in 
each of these. The first sweep at 1410 EST shows a 
fairly regular scatter trace which originates at the 
second-order F2 vertical-incidence trace and exhibits 
the separation discussed above. There is another fine 
line of scatter echoes coming off the second-order F1 
vertical-incidence trace which has been demonstrated 
to be ground scatter propagated by the F1 layer. In 
the next sweep at 1458 the F1 layer is starting to dis- 
appear and the scatter propagated by it is weaker. In 
the last sweep at 1652 the F1 layer and the scatter asso- 
ciated with it have disappeared,—also the apparent 
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skip ranges for F2-propagated scatter are noticeably 
longer than in the previous sweeps, because of the late- 
ness of the hour. 

A method of determining the source of scatter propa- 
gated back to the receiver via any layer was first shown 
by Dieminger.?:4 The simple proof of the technique ap- 
pears in a paper by Peterson.® If the ordinary wave 
sweep is plotted on linear co-ordinates, a line from the 
origin (zero frequency and zero height) tangent to the 
second-order trace marks the locus along which the 
minimum-time-delay backscatter echoes of the ordinary 
wave will fall as the frequency increases, if these echoes 
are from the ground. If they are from the distant E 
layer they should lie along a line starting at a point 
about the E-layer virtual height below the tangent 
point. The technique is correct for essentially transverse 
propagation. 


1400 


1200 


1000 


800 


600 


SLANT RANGE, KM 


400 


1 


TIME DELAY 
TRACE 


POINTS OF MINIMUM 


** ON BACKSCATTER 
aE 


200 


8 
MEGACYCLES 


(A) TANGENT-LINE THEORETICAL LOCUS OF LEADING EDGE OF Fo- PROPAGATED 
GROUND SCATTER ECHOES 

(B) SIMILAR LINE FOR F; LAYER 

(C) THEORETICAL LOCUS FOR Fp -PROPAGATED ECHOES FROM DISTANT E LAYER 


(D) THEORETICAL LOCUS FOR Fp- PROPAGATED GROUND -SCATTER ECHOES AT 
THE SKIP RANGE 


Fig. 4—Analysis of multifrequency backscatter echoes following sec- 
ond-order F1 and F2 reflections at Sterling, Virginia, January 23, 
1953, 1124 EST. 


Fig. 4, for January 23, 1953 at 1154 EST which is 
similar in principle to figures shown by Dieminger and 
Peterson, is presented as an illustration of the tech- 
nique. The solid curved lines are a plot of the second- 
order F1 and F2 layer ordinary-wave traces, and the 
crosses are a plot of the leading edge of the backscatter 
echoes. A straight line (A) drawn from the origin tangent 
to the F2 ordinary-wave second-order trace falls close to 
the line of the backscatter echoes at the outset. The line 
of echoes represented by the three lowest crosses illus- 
trate the separation often seen and appears, as men- 
tioned above, to be associated with the extraordinary 
wave. The crosses at 7.6 and 8.0 mc represent a jog 
often seen in this type of trace and could represent a 


5 A. M. Peterson, “The mechanism of F-layer propagated back- 
scatter echoes,” Jour. Geophys. Res., vol. 56, pp. 221-237 ine! 1951. 
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region at which resolution of the two separated echo 
lines ceases to exist. Similar but larger jogs sometimes 
seen on disturbed days are likely to be due to steep 
horizontal gradients in the ionosphere. 

Another straight line (B) from the origin tangent to 
the F1 layer ordinary-wave second-order trace falls 
very close to the points of the type identified with F1 
propagated ground scatter in the discussion of Fig. 3. 

The isolated slant line (C) below the F2 tangent line 
is the line along which the F2-layer-propagated scatter 
echoes would be expected to fall if they were scattered 
back from an E layer at a virtual height of 105 km 
above the ground, and it certainly appears here that 
the scatter does not come from the £ layer, especially 
at the start. The line of scatter echoes approaches 
nearer to the E-layer line at greater distances, but here 
it must be remembered that ionization at the reflection 
point is not the same as overhead ionization and the 
more distant scatter could have been propagated by a 
region of greater ion density. 

Peterson® introduced the concept of “time-delay 
focus” and demonstrated that it causes the amplitude 
of a scatter echo to rise rapidly after the minimum-time- 
delay echo, which represents waves scattered from a 
greater distance than the skip distance, to a maximum 
at a delay time still short of the skip-distance delay 
time. In Fig. 4 the dashed line (D) above the minimum- 
time-delay line marks the delay time for the skip- 
distance ray, calculated by the use of N. Smith’s trans- 


_mission curves for the curved earth and curved iono- 


sphere.® It will be noted that the difference between the 
minimum-time-delay ray and the skip ray begins to 
become negligible at about 800-km slant range. Beyond 
950 km the separation apparently widens again. How- 
ever this is the case only because the minimum-time- 
delay tangent line was drawn*on the flat-earth, flat- 
ionosphere hypothesis and actually should curve up- 
ward and approach the skip ray asymptotically. The 
upward curvature is caused by the fact that a given 
apex half angle, angle ¢o on the ray trajectory, is asso- 
ciated with a greater slant range in the curved-earth 
case than in the flat-earth case, the relative increase 
becoming greater as distance increases. 

Out of a dozen sweeps for the period under con- 
sideration which were analyzed in this way, only one 
yielded a disposition of points which might be con- 
sidered to be in disagreement with the proposition that 
the short-range F2-propagated scatter comes from the 
ground. 

Considering the disturbed period, Fig. 5 shows a 
sequence of eight records for December 17 beginning 
just after local noon. The records are characterized by 
high absorption of the lower frequency echoes and 
great distortion of the traces with distant echo traces 


_ °N. Smith, “The relation of radio sky-wave transmission to 
Bice yoke measurements,” Proc, I.R.E., vol. 27, pp. 332-347; May, 
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splitting in odd ways, the entire configuration of the 
sweep changing rapidly from sweep to sweep. Nearly 
horizontal traces with sharp jogs and straight-line in- 
clined traces in a regular pattern are due to pulse 
interference. 
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Fig. 5—Sweeps for December 17, 1952. 


Although few vertical-incidence sporadic-E reflections 
are seen on this set of records, many were seen on 
vertical-incidence records made at the same time at 
nearby Ft. Belvoir, Virginia. The Sterling records for 
the following day, December 18, showed a great deal 
of sporadic-E reflection. These records will be discussed 
-below under sporadic-E phenomena. 

In an analysis of a dozen of the records for this period 
by the tangent-line method, six left no doubt but that 
the backscatter near vertical incidence was from the 
ground. The other six were somewhat doubtful but only 
one of the six gave a distinct impression that the scatter 
might be coming from the £ region. 

To test the merit of using such backscatter records 
for determination of the maximum usable frequency 
over a fixed distance and in particular to test the 
validity of the assumption that backscatter takes place 
at the ground rather than at the £ layer or other 
ionospheric region, use was made of vertical incidence 
records obtained during these same two periods at the 
ionosphere station at Batavia, Ohio, the midpoint of 
the 1,150-km path between Sterling and St. Louis. As 
mentioned previously the backscatter records were ob- 
tained using rhombic antennas pointed toward St. 
Louis. 

Unfortunately the transmitting operations near St. 


- Louis had been terminated when data-taking was begun 


in this experiment; however, close agreement between 


the observed delay times and muf for oblique trans- 
mission over the path and the delay times and muf 


derived from the midpoint records had already been 
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well established. Vertical-incidence sweeps made at 
Batavia every three minutes were scaled by use of 1,150- 
km transmission curves to obtain the muf and the 
virtual height. Using these virtual heights, path-delay 
times were obtained from geometrical calculations 
under the assumption that the scatter was from the 
ground, and then the apparent muf was scaled at these 
delay times on the 12-minute backscatter sweeps, ap- 
propriate sweep-time corrections being made. In the 
few cases where the scatter records had two branches 
at the correct delay time, scalings from both were 
plotted. However, in the subsequent analysis of results 
only the lower branch was considered, unless the upper 
branch showed evidence of being a much stronger echo, 
in which case it alone was considered. 

Use of the geometric method, to calculate the delay 
time for the equivalent triangular path, was justified 
by previous results of analyses of round-trip travel 
time for the 1,150-km path which indicated that the 
geometric method gave travel times agreeing with ex- 
periment within about 2 per cent, or almost within ex- 
perimental error. It should be noted that our “slant 
range” is really a virtual slant range, since it applies to 
the equivalent triangular path. 
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Fig. 6—January 22, 1953; 1,150-km path, assuming ground scatter. 


Fig. 6 is for one day of the quiet period, January 22, 
1953, dots indicating the path muf computed from the 
midpoint data and crosses indicating the path muf 
computed from the backscatter. It is evident that the 
two sets of muf usually agree to within a few per cent 
and there is no doubt but that the scatter originates at 
the ground. The small systematic difference is con- 
sistent with results obtained in the oblique incidence 
experiment between Sterling and St. Louis mentioned 
previously. As a further check, the small rectangle at 
1726 EST indicates the failure of 10-mc transmis- 
sions from Station WWYV, at Beltsville, Maryland, ob- 
served at St. Louis. It is seen that this falls right on the 
graph of the scatter data. A similar result was obtained 


on January 23. 
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Fig. 7 is for one day of the disturbed period Decem- 
ber 18, 1952, the dots and crosses having the same 
meaning as before. The results here are vastly different 
from those for the previous period. The mut derived 
from the backscatter sweeps, assuming the ground to 
be the origin, have much larger values than those de- 
rived from the midpoint data. 
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Fig. 7—December 18, 1952; 1,150-km path, assuming ground scatter. 


When these results were noted the next logical step 
appeared to be to assume that the distant E layer was 
the source of backscatter. A height of 105 km was as- 
sumed for the E layer and to simplify calculations, a 
constant F2-layer virtual height of 300 km was used. 
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Fig. 8—December 18, 1952; 1,150-km path, assuming E-layer scatter. 


The amount by which all slant ranges calculated for 
ground scatter should be diminished for this case is 
approximately 240 km. Typical results are shown in 
Fig. 8. Here it is evident that although there is consider- 
able variation in the difference between the midpoint 
and backscatter results, the agreement on the whole is 
very close. 

Frequency diagrams of number of observations vs 
per cent difference (backscatter frequency minus mid- 
point frequency), (a) for the January period assuming 


ground scatter, (b) for the December period assuming 
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ground scatter and (c¢) for the December period assum- 
ing E scatter, appear in Fig. 9. Results are expressed in 
‘Tablet: 


TABLE I 
b Per cent Difference 
Assumed er cS Median Exceeded in 
Period Scatter Oneerras Difference 
A Per cent 95% 5% 
pare Boe of cases of cases 
Tanuary 21-33, 
3 ground 73 4.4 —0.8 11 
December 17-19, 
2 ground 73 22.2 13e7 3959 
December 17-19, 
1952 E layer 95 1.8 —3.2 21.2 


The median value of 4.4 per cent by which the values 
of muf for January 21-23, 1953 determined from the 
backscatter, assuming ground scatter, exceed those de- 
termined from the midpoint data, should be compared 
with results of the oblique-incidence, two-way pulse 
experiment previously conducted over the same path 
in which the median value of the per cent difference be- 
tween the muf of the actual transmitted pulses and the 
muf derived from the midpoint data for daytime in 
January 1952 was 3 per cent. Thus it would appear that 
for the conditions of January 21-23, 1953 the backscat- 
ter assumed to be from the ground was an indicator of 
muf to within 1 or 2 per cent error. 

The median value of 22.9 per cent by which the 
values of muf for December 17-19, 1952 determined 
from backscatter, assuming ground scatter, exceed 
those determined from the midpoint data makes the 
ground-scatter hypothesis seem invalid. However, when 
the backscatter muf is based upon assumed scattering 
from an E layer about 105 km high, a median difference 
of only 1.8 per cent results, suggesting that scatter by 
the distant E layer was the source of most of the echoes 
and that perhaps the correction should be smaller, as 
would be the case for a slightly lower height of the E 
layer. 

Also, in Fig. 9 the similarity between the distribu- 
tions for the January 21-23 period, assuming ground 
scatter, and the December 17-19 period, assuming E 
scatter, is noted. The wide spread of differences of a 
small number of points in the December period is prob- 
ably due to the fact that this period was characterized 
by stratifications and turbulence. Violent fluctuations — 
of F2 layer muf with time like those on December 17. 
indicate space fluctuations of a similar magnitude, 
which would enable the relatively broad-beamed an- 
tennas used in this experiment to find paths of equiva- 
lent delay time having substantially higher muf than 
that-over Batavia. ‘ 

One is led to believe by the comparison with the 
Batavia records that the scatter in the December 
period was returned from the distant E layer,—yet 
the scatter trace at lower frequencies merges with the 
second-order F2 vertical incidence reflection indicating 
that short distance scatter was from the ground! It is. 
conceivable that there could have been an E layer struc- 
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ture from which the scattering was much more promi- 
nent at oblique incidence than at vertical incidence; but 
if this were the case, the transition from ground scatter 
to E-layer scatter should appear on the December 17—19 
records, both being visible on some part of a sweep. No 
clear-cut transition is visible on any of the records. 
Jogs in the scatter trace like that of Fig. 4 at 7.8 mc 
could be construed to indicate a mode change, but the 
effect of F2-layer turbulence may also have played a 
dominant role. The latter, by producing irregular slant 
range records could probably mask any transition 
effects. 
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Fig. 9—Distribution of muf differences. 


DIRECT OBLIQUE REFLECTIONS FROM SPORADIC E AND 
GROUND SCATTER PROPAGATED VIA SPORADIC E 


The sweep-frequency records tend to bear out con- 
clusions previously reached in studies of range-time 
records at a fixed frequency that certain echoes repre- 
sent direct, oblique reflections from sporadic-E regions 
while others represent ground scatter propagated by 
sporadic E£.’ 

The first of two sweeps in Fig. 10 for December 19, 
1952 at 0450 was chosen for examination, first, because 
at that hour the F2-propagated ground scatter is not 
seen and, second, because the F2-layer first echo is in a 
position where it is not mixed with the sporadic-E 2nd 
echo and the ground scatter propagated by sporadic E. 
In this sweep the first F2 ordinary-wave echo can be 
seen faintly in the lower left-hand corner with a mini- 
mum virtual height of about 320 km and critical fre- 
quency of 3.2 mc. Below it at a little over 100 km is 
seen the first-order sporadic E reflection terminating at 
about 4.2 mc. 

Above the regular sporadic E reflection is clearly seen 
the diffuse second-order reflection. Out of this is seen 


7 W. L. Hartsfield and R. Silberstein, “A comparison of cw field 
intensity and backscatter delay,” Proc. I.R.E., vol. 40, pp. 1700- 
1076; ee at 1952. 
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to develop a very dense echo which recedes to 800 km 
at 19 me. A fainter echo group finally fades out at 21 
mc and 1,150 km, The strong echo appears to be ground 
scatter propagated via sporadic E ionization since it 
grows out of the second-order echo. It is also evident 
that, since the sporadic-E distance factor for an 800-km 
slant range is about 2.8, the maximum frequency value 
at this distance would have to be associated with a 
sporadic £ region at a control point where the maximum 
frequency value, assuming power dependence is not 
important, is of the order of 7 mc. Below the dense 
echoes there are weak echoes running out to about 13 
mc at 300 km which could be due to oblique reflections 
directly from a localized sporadic E region. 

In the sweep for 0514 on the same day the F2 reflec- 
tions are not seen. Here the apparent direct reflections 
from sporadic E go out beyond 18 mc at about 400 km. 
The dense, longer-range echo again appears to be 
ground scatter propagated via sporadic E. 

An oblique echo which exhibits a fixed range with in- 
creasing frequency such as the sharp one at 1,150 km 
in the sweep for 0450 probably represents reflection or 
scattering off a small dense region of sporadic E or 
ground scatter propagated by such a region. 


Fig. 10—Sweeps for December 19, 1952. 


Fig. 11 is a series of three daytime sweeps made on 
December 18, 1952, which illustrates phenomena asso- 
ciated with sporadic E in addition to the regular scatter 
traces. In the sweep for 1155 EST direct sporadic E 
scattering from an obliquely situated region is evident 
out to nearly 11 mc at 200 km. Another heavy sporadic 
E trace appears to grow out of the second-order spo- 
radic-E echo and blots out the first-order F2 reflections, 
finally terminating at 13 or 14 mc and 500 km. This 
appears to be the ground-scatter type because of its re- 
ceding range and association with second-order spo- 
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radic E. The sweep for 1231 shows a similar effect appar- 
ently with direct reflections from many small clouds and 
stratifications of sporadic E at short ranges. 

In the sweep for 1307 the scatter appearing after the 
first-order F2 reflection upon casual inspection might be 
interpreted as direct scatter from the F2 layer since the 
dense part starts after the F2 ordinary-wave critical 
frequency, continuing to 9 mc, and since a similar ap- 
pearing and equally dense trace occurs at about twice 
the delay time, apparently superimposed on the regular 
scatter trace emanating from the second-order F2 trace. 
A few other sweeps around this time had similar fea- 
tures, but there seems to be no definite reason why direct 
F2 scatter should have appeared. 

Questions arising out of study of the record for 1307 
led to a determination of how a particular mode, the 
“1F” mode proposed by Eckersley® would appear. This 
mode is the case of the wave being reflected at the F2 
region, striking a distant scattering source in the E layer 
and returning directly to the transmitter. The reverse 
passage of the wave would be equally possible. Results 
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Fig. 11—Sweeps for December 18, 1952. 


of an analysis indicate that the mode grows out of the 
first-order F2 reflection a little below the critical fre- 
quency. It should go to higher virtual heights than the 
heavy scatter emanating from the first F2 echo in any 
of the sweeps in Fig. 11. The thin steep trace going from 
about 6.6 mc at 470 km to 7.8 mc at 600 km in the 
record for 1307 was thought possibly to be such a trace, 
but does not agree well with the calculated curve. Be- 

8 T. L. Eckersley, “Analysis of the effect of scattering in radio 


transmissions,” Jour. IEE, vol. 86, pp. 548-563, disc. 563-567 
Fig. 23b, p. 559; June, 1940. : 


cause of the low angle, the phenomenon, if seen, should 
be chiefly for fairly short distances, and it might be 
deduced from the records so far examined that the mode 
does not play a very important role. 
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Fig. 12—Sweeps for December 24, 1952. 


Fig. 12 is a sequence of fourssweeps for December 24 
with a 1,000-km range scale illustrating more sporadic- 
E phenomena in addition to the regular scatter growing 
out of the second-order F2 modes. Echoes appear which 
have the properties of direct and oblique reflections from 
sporadic-E regions and ground scatter propagated by 
sporadic E in rapidly changing configurations. “M” re- 
flections are also seen. 


CONCLUSION 


Sweep-frequency backscatter records, by introducing 
the frequency parameter into scatter-sounding tech- 
niques, make it possible to obtain a large number of 
samples for study of phenomena on a given frequency 
or at a given distance in any one day. Also the develop- 
ment of echoes as frequency increases gives strong clues 
as to their origin and as to the relative quietness or dis- 
turbance of the ionosphere in the period under study. 
Under certain quiet ionospheric conditions the meas- 
urement of maximum usable frequency and therefore 
of skip distance, from backscatter records under the as- 
sumption that scattering is from the ground yields ex- 
cellent results; however, it was demonstrated that the 
same assumption yields very erroneous results during a 
period of intense sporadic E reflection. Under these lat- 
ter conditions the distant scatter appears to originate 
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at the E layer, although nearby scatter still appears to 
originate at the ground. To date no clear transition be- 
tween the two sources has been observed on any one 
sweep and not enough data have been taken to afford 
observations on the change between the type of iono- 
spheric conditions where distant scatter appears to come 
from the E layer and the type where it appears to come 
from the ground. 

Relative to the effectiveness of backscatter measure- 
ments in skip distance indication it can be said that a 
single frequency technique may fail during a low- 
latitude disturbance, which is not forecast at present 
by forecasting agencies; yet during such disturbances 
signal behavior is erratic and skip-distance information 
would be desirable. It is noteworthy? that commercial 
equipment for skip-distance determination permitting 


§L. C. Edwards, COZI—Communications Zone Indicator, Elec- 
tronics vol. 28, pp. 152-155; August, 1953. 
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measurements on several frequencies has recently been 
described in the literature. 

The sweep-frequency technique tends to confirm the 
presence of direct scatter from discrete regions (“clouds”) 
producing sporadic-E reflections and of ground scatter 
propagated by sporadic E but the echo patterns are 
frequently so complex as to make a positive statement 
about any one very difficult. 

The sensitivity of the records to varying conditions 
of the ionosphere indicates that disturbance tracking 
and forecasting by its means may be possible.7:! 

The use of high power with very high-gain, narrow- 
beam antennas should make a study of night F2-layer 
propagated backscatter echoes possible, thus affording 
information on propagation of the lower, high fre- 
quencies during the hours of darkness. 


10 R. S. Lawrence, “Continental maps of four ionospheric dis- 
turbances,” Jour. Geophys. Res., vol. 58, p. 219; 1953. 


Some Stochastic Problems in Wave Propagation—Part II* 


JOSEPH FEINSTEINT 


Summary—tIn part I some of the properties of wave energy re- 
flected from surfaces subject to random height variations were in- 
vestigated in the previous issue. 

Here we ascertain the effects of refractive index fluctuations 
within a volume upon the properties of waves traversing the medium. 
The results are applied to problems encountered in tropospheric and 
ionospheric wave propagation. 


VOLUME PROBLEMS 
Tropospheric Applications 
R uecrea WORK ‘' on the effects of refractive index 


fluctuations upon wave propagation has been 

based upon the primary scattering produced by 
the dipoles equivalent to the medium fluctuations. An 
alternative viewpoint, developed here, considers the 
phase incoherence introduced into the wave fronts by 
these refractive index fluctuations. This approach takes 
account of multiple scattering which, although it may 
play an important role in the ionosphere is not likely to 
be of significance in the troposphere. While the two 
methods are equivalent when only single scattering 
enters, the phase incoherence approach leads in a natu- 
ral way to a diffraction phenomenon which has thus 
far been overlooked in the scattering theory treatment 
of propagation beyond the horizon. In addition, the 
determination of the fluctuation characteristics of wave 
energy within line of sight seems better adapted to the 


methods of this paper. 


Phase Fluctuations Produced by a Finite T. hickness 
Screen: We consider the field produced at plane Q, a 


* Revised manuscript received by PGAP, January 18, 1954. 

+ Bell Telephone Labs., Murray Hill, N. J. : ’ 

1H. G. Booker and W. E. Gordon, “A theory of radio scattering 
in the troposphere,” Proc. I.R.E., vol. 38, p. 401; April, 1950. 
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Fig. 1—Single slab geometry. 


distance D, from a slab of medium of thickness R, il- 
luminated by a transmitter at a distance D, (Fig. 1). 
Employing usual second-order phase approximation, 
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27D,D>2 —o 
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for R&Dj 2 where Xo, yo are co-ordinates in the xy plane 
of Fig. 1, the origin being taken as the intersection with 
the xy plane of a perpendicular from the transmitter; 
Ad is the phase departure from its mean value intro- 
duced by the fluctuations in the slab. The association of 
this phase with a single point in a plane carries with it 
the implication of phase integration through the slab, 
Ao=kfrAu(n, t)drn, each trajectory being specified by 
its point of intersection with the reference plane. For 
simplicity, Au will be taken as real, thus neglecting ab- 
sorption. We are interested in obtaining the correlation 
pattern at Q, defined as: 


Co(x0, Yo; 5, v) 


1 T 
= lim —f[ dt- Eg(xo, yo, t) Eq*(xo + s, yo + v, t). (2) 
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Using (1), and interchanging order of integration: 
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To proceed further, the statistical properties of Ad 
must be specified,—in particular the probability dis- 
tribution of “Ad(x, y, t) -Ad(x’, y’, t).” Since each Ad 
is composed of a large number of random contributions 
as a result of the integration of (Au) over each trajec- 
tory, we may invoke the central limit theorem to justify 
the assumption of a bivariate normal joint distribution. 
Then the difference of the Ad’s is also normally dis- 
tributed with variance (A¢)?(x—x’, y—y’), a function 
only of the difference of the coordinates if the fluctua- 
tions are homogeneous. Average over time may now be 
replaced by average over the probability distribution: 


1 oy 
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= e—3(A¢d)*(x—2', y—y’), (4) 
Making the transformations £=x—x’, »=y—y’ in (3), 
and utilizing the Fourier integral theorem: 
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For a plane wave (Dix) at normal incidence 
(x9 =yo=0) we continue to obtain the field correlation 
pattern at the slab, over all parallel planes. This result 
has been obtained previously for a one-dimensional 
screen.2 In the more general case, a scale factor enters, 
and an oscillating phase factor is superimposed upon the 
pattern. We next determine the mean square difference 
of phase along two trajectories, 1 and 2. 
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If we define a correlation function in the medium by: 


(6a) 


1 Tr 2 ee 
aie BPG t)-Au(r’, t) = (Au)?p(7, 7’) (7) 


then for parallel trajectories with perpendicular spacing 
“q,” (6a) becomes: 
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where the transformation s=7,;—7;'=fr2.—7re' has been 
made; for isotropic turbulence, p(s, ¢) =p(W/s?-+q). If 
the slab thickness, R; is several times the scale of turbu- 
lence, limits of + © may be employed on the integrals 
involving p, and (R—s)<R. Finally then, for a correla- 
tion function of the form p(u) =e“, one obtains: 


(Ad), = 2k(An)™ R-l-/a(1 — ew") (9a) 
while with p(w) =e7!"l/! 
(Ad)? = 4k?(Ay)?- R:l- E — _ (21) (9b) 


where K, is the Bessel function of the second kind for 
imaginary argument. For substitution in (4) and (5) one 
of course sets g= V/(x—x’)?+(y—y’)?. 

The technique of phase integration is essentially a 
geometric optical one. Consequently diffraction effects 
are neglected over the thickness of the slab, which is re- 


. © F HIG. Booker, J.A. Ratcliffe, D. H. Shinn, “Diffraction from an 
irpeler screen,” Phil. Trans. Roy. Soc. A, vol. 242, p. 579; Sept., 
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quired to have a value a few times the scale of coherence 
in virtue of the assumptions above. This is analogous to 
neglecting diffraction over a distance of a few slit 
widths in front of a slit, an approximation which does 
not appear serious. 

Extension to multiple screens: To consider the effect of 
a series of slabs upon the wave, we take the field in- 
cident upon slab n at distance z, from the source as that 
which has emerged from slab (7 —1) at 2,1, and evalu- 
ate the additional phase incoherence introduced on a 
plane at 2n41. 
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The correlation pattern at this position is then: 
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After transforming variables through §=x—x’, n=y 
—y’, the time average of the product of the fields at 2, 
may be replaced by the correlation pattern (5), after 
setting x»=x, yo=y, S=E, v=, Di=2n1, D2=2%n—-2na 
in this pattern. Thgtime average over the phase fluctu- 
ations introduced by the mth screen is taken separately, 
since events on this screen are not correlated with the 
incident wave. This yields an expression analogous to 
(4). The x dependence in (11) has the form 
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so that integration over x yields a delta function of 
argument [(en+1/zn)E—s]. Integration over & then 
amounts to setting £=sz,/(Zn41). (This procedure is 
identical with employing the Fourier integral theorem, 
as was done in deriving (5)). Analogous results follow 
for y and 7, so that (11) yields: 
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For a group of m discrete slabs positioned at distances 
Zm from the source, the correlation pattern at distance D, 
is given by: 

ene 
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This result is subject to a rather simple interpretation. 
Since the product of a group of exponentials is equiva- 
lent to a sum of exponents the phase incoherence effects 
introduced by the slabs, each modified by a scale factor 
in accordance with the ratio of source-slab to source- 
receiver distance, add at the receiver plane. Precisely 
this result would follow from the geometric optical 
method of integrating the phase along the total path of 
each ray from source to receiver. There is this differ- 
ence, however. The solution obtained here by the meth- 
ods of physical optics is a consequence of the infinite 
limits on the various integrals taken over the x —y plane, 
i.e., each slab of medium was taken to be infinite in 
the transverse directions. Geometric optics does not im- 
pose this requirement. So long as diffraction effects are 
unimportant over the entire path length a one-to-one 
correspondence exists between points on the wave-front 
at difference z, rendering the concepts of geometric 
optics valid. This is the situation for 1/A— ©. Once dif- 
fraction enters however, it is an area on a given wave 
front which contributes to the field at a point on a suc- 
ceeding wave front. The significance of (13) lies in its 
statement that the scattering of energy out of the direct 
ray under these conditions is exactly compensated by 
the diffracted energy from adjacent portions of the wave 
front. It should be noted that if the angular range over 
which significant diffraction occurs becomes sufficiently 
large to make it necessary to take account of the varia- 
tion of the distance attenuation and obliquity factors 
assumed constant in (J), these results no longer hold. 
For the usual values of “1” encountered in the tropo- 
sphere this condition appears unlikely. 

We should now like to generalize (13) to continuous 
media. To this end we may visualize the distance D 
divided into sections R, each at least several times the 
scale of turbulence, J, so that the coherence between 
sections is negligible compared to that within a section. 
Then provided D>>R so that the scale factor remains 
essentially constant within a section and neglecting for 
the moment the apparent failure of the second order 
phase approximation as the distance between slabs be- 
comes small, we may substitute an integration of the 
exponent in (13) for the discrete set of product factors. 
Neglecting phase factor (we may take xo=yo=0; the 
Fresnel zone argument k(s?+v*)/2D is generally slowly 
varying compared to correlation pattern amplitude), we 
obtain for gaussian correlation function: 
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Co(s, v) = (14a) 
while utilizing (9b), the exponential correlation function 
yields: 
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For homogeneous turbulence, (Ay)? and J may be re- 
moved from the integral, so that (14a) reduces to: 
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We must now consider the significance of these ex- 
pressions in the light of the fact that the second order 
phase approximation upon which they are based is 
valid only for energy transmission from section to sec- 
tion within cones of small angle; except for cases wherein 
geometric optics is applicable, this condition would 
appear to be violated as the continuous limit is ap- 
proached. A quantitative error criterion may be obtained 
from a power series expansion of the exponential in 
(15). The terms in such an expansion may be interpreted 
as the various orders of scattering; the contribution of 
single scattering, for example, is given by the linear 
term and represents the sum of the scattered waves 
produced by all the slabs, the effect of each being calcu- 
lated in the absence of all the others. It is apparent 
that our phase approximation is valid for this term, ex- 
cept possibly for a few slabs at each end of the path. 
Similarly the quadratic term gives the second order 
interactions which number n(n—1)/2 for n slabs, of 
which only ” are between adjacent slabs. For D>, 
consequently, the approximation will break down only 
at quite high orders of scattering, and under conditions 
such that these high orders are of importance, several 
of our other approximations will fail. For the usual con- 
ditions encountered in tropospheric propagation only 
single scattering is of importance so that the expansions 
of (15a, b) need be carried only through the linear 
terms. 

Information may be obtained about the frequency 
power spectrum of the received signal, by forming the 
temporal (auto-) correlation: 
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To proceed through steps analogous to (6), (7), and 
(8) we must define a new type of correlation function: 
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If it is assumed that “Ad(x, y, t) and Ad(x’, y’, +7)” 
continue to have a bivariate normal distribution then 
(4) is applicable. Finally we must make some assump- 
tion regarding the space-time correlation function in 
(17). If we write 
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then (16a) becomes: 
P(r) ~~ e-al—p(r)]. giwor (19) 


where a=bk?(Au)7D. “b” having the value 7m for a 
gaussian “p,” and 2 for a brownian “p.” 

The power spectrum is now obtained as the transform 
of the autocorrelation function: 


P(w) = fi Po(s)-dr. (20) 


For p(r) =e-7/27:, we obtain 
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Diffraction Beyond the Horizon Produced by Phase In- 

coherence: In this section we consider the effect of phase 

incoherence due to turbulence in modifying the dif- 

fraction field which results from earth curvature. Fig. 2 
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Fig. 2—Diffraction beyond the horizon. 


illustrates the physical situation. The new element 
here is the fact that the point of stationary phase in 
the x-co-ordinate, x =0, no longer lies in the volume of 
integration. We obtain the field at Q by summing the 
contributions which arise from the field on the mid- 


plane M, each with its proper phase. 
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where normalization with respect to the free space field 
has been introduced; A@ is the phase variation intro- 
duced into the wave in going from M to Q. The average 
power received at Q is then: 
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After usual change of variable £=x—x’, n=y—y’, 
correlation pattern given by (13) with amplitude as in 
(14b) is substituted for EyEy*; argument (x0, yo; s, v) 
of the latter is of course replaced by (x, y; &, 7). For the 
time average of the difference of the (A@)’s is (22) we 
gain employ (4), (9b), and finally (14b): 
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The integration over y yields (tD/k)/5(n), so that 
integration over 7 merely consists of setting n=0. 
The factor F?(x) has been introduced to connote the 
fact that between x=a, the surface of the earth, and 
x =X», the horizon, the field increases gradually from a 
very small value to its free space value. Since 4k?(Ay)*/D 
<1 for the usual tropospheric case, we may employ the 
first two terms of the exponential expansion: 
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where homogeneous turbulence has been assumed. For 
(Au)?=0 over the whole path the integration of (24) 
would give rise to the small diffracted field at Q char- 
acteristic of propagation over a spherical earth. Since we 
are interested in the contribution which results from 
phase incoherence we shall neglect this first term and 
begin the X-integration from x9. In addition we may 
neglect the phase variation produced by exp (—tké?/D) 
compared to exp (2zkx£/D) provided the vertical dimen- 
sion (x) of the effective contributing volume is much 
greater than “J,” the scale of coherence (which de- 
termines the effective limits of £). With these changes 
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where g=2kixo/D. For g>>1, corresponding to large 
scale turbulence, there results: 
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where R, is the radius of the earth. It is interesting to 
compare this value with the contribution made by di- 
rect scatter from turbulence lying within line of sight. 
For the exponential type correlation function used 
above, an expression for the average power per unit 
volume has been given in reference (6). For large “g,” 
integration of this expression over the trough shaped 
region above line of sight yields precisely the same de- 
pendence upon the parameters as does (26); the mag- 
nitude is virtually identical with that of (26) for the 
same (Ay)?//. In practice furthermore one would expect 
this turbulence parameter to decrease with increasing 
altitude. Since (23) indicates that it is a weighted mean 
of this parameter over the total path which is to be 
employed in (26), whereas only its value above line of 
sight is to be used in evaluating the direct scatter, the 
contribution of (26) can be expected to strongly dom- 
inate over the direct scatter at large distances. From the 
point of view which treats the effects of medium fluctua- 
tions as scattering, the field contribution derived here 
represents the diffraction of forward scatter. 

To obtain correlation pattern over a vertical plane at 
the receiver, we add a displacement (s, uv) to (x, y) in 
phase term of (22). Integration over y then leads to 
5[n—(v/2)], so that we set »=v/2 to integrate over 7. 
Integration over & can be performed next, leaving: 
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This expression can be evaluated for g>>1. 
For the special cases where s or v equal zero, we ob- 
tain for the normalized correlation amplitudes: 
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where ap=x0/D. Corresponding expressions for direct 
scatter are similar, although not identical, but the differ- 
ences do not seem sufficient to admit of an experimental 
discrimination. What experimental evidence there is 
seems to be in agreement with forms such as (28). 

Too, these correlation expressions prove of value when 
effect of objects adjacent to source or receiver can be 
represented by means of images. Example: If effect of 
ground reflection upon transmitter and receiver antenna 
fields can be represented by images at positions separated 
in height by hr and he respectively (distances above 
ground of 447 and hp), results of (28) and similar devel- 
opments lead to formula below for average power: 


P(hr, hr) — Pol4 — 4Cy(hr) cos khe — 4Cy(hr) cos khr 
+ 2Cy(hr == hr) cos k(hr + hr) 
+ 2Cy(hr + he) cos k(hr — he)). 


where Pg is given by (26) and Cy is the normalized 
correlation amplitude in the vertical direction. 


Tonospheric Applications 


This class of problems is characterized by the de- 
pendence of the degree of inhomogeneity upon the 
mean value of refractive index when the state of 
turbulence is specified in terms of electron density 
fluctuations, the physically significant variable. In addi- 
tion, it differs from problems of the type considered in 
section (a) because of the large variation in mean re- 
fractive index encountered by a wave, leading generally 
to total reflection at some level. 

Formulation of Ionospheric Fluctuations: The phase 
integral representation of a field reflected from an in- 
homogeneous medium is given by? 


tkxo sin 6+ 22k 
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where x) is the total horizontal distance between 
transmitter and receiver, u is the refractive index with 
mean value a function of height z, u(zm) =sin 0, and @ is 
the angle made with the normal by a given trajectory 
at take-off. I 

We set w=o+Au, and expand the phase to first 
order in (Ay): 
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The method of Sraconery phase is now applied to 
the mean terms, to find the trajectory angle 49, which 
corresponds to a given range x». 
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Utilizing the second order approximation to the phase, 


(29) can now be written: 
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The average power takes the form: 
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We now make our usual assumption that the differ- 
ence of the two random phase fluctuations is normally 
distributed. Some justification for this assumption exists 
in this case regardless of the distribution of (Au) at 
any given point, in virtue of the central limit theorem, 
since each integrated random phase is the resultant of a 
large number of independent contributions. We proceed 
to determine the variance of this phase difference: 
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It is now desirable to introduce the physical variable, 
electron density, related to » through: 
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and replacing linear distance in the argument of the 
correlation function by As, arc length along the trajec- 
tory, an approximation which is justified provided the 
radius of curvature of the trajectory is large compared 
to the scale of coherence, (35) becomes: 
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where the transformation ¢=z—2' has been effected, 
and z has been set equal to z’ except in the correlation 
function; € is the lateral displacement of trajectories 
specified by @ and 6’ at height z, so that the argument 
of the second correlation function is the distance be- 
tween points at height z on trajectory 0, and 2z’ on 6’. 
We now determine (As) and ¢ is terms of z 
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» 9%0(z) 
E = x0(z) — xe(z) = (6 — 6’) org Se BF(6, z) (42) 


where we have set B=0—6', F=0x/00. 
Returning to (38), for a correlation function of the 
form p(u) el the integral over ¢ becomes: 


‘ ee it £10” 2£¢ sin 0 | 
f aye P-(uc?—sin? 6) acs 6) a a5 ny 2 sin? a/ po? —sin? 6 


; ae 2 sin? 6 
V Mo” — Sin {1 he rete, ; ) 


Ho 


=/rl (43) 


where infinite limits have been employed, since the 
thickness of the ionosphere is much greater than the 
scale of coherence. 


R2ly\/ 3 1}? 


(~) f te (0? 
2 N 0 wo? V wo? — sin? 6 


(Ag)? = 


_ BF? (6, z) ak sin? *) 
{1 —¢e ? p Sods (44) 
The average power takes the form: 
i" (80) 
(268 — 2608 — 6”) ——e 
Ps, al do-dB-e ? ; -e—2(Ad)p?. (45) 


To obtain the correlation pattern on the ground, we 
add a displacement s to x». This has the effect of alter- 
ing 0) by an amount, Ad» given by: 


lasing, [=| (65) 
anne | Tass 


which we may write as s=(A6o)-g(90), where g(40) is 
the derivative on the right hand side. The change in 
$''(00) is of second order, so that the phase term in (34) 


takes the form : 
Ss 2 
[v-we-(r-n2)] 
en, * 


s = Ax = 


'’(80) 
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This gives rise to additional terms 
Os iL / B\e 5 Gas 
6") + ( \+ | 
g 2\g i 
Then 
P,,(s) -{f dé. dg 
i Ris as le fas. ie 1 
at a) [o(8 =) AG a) ale) 78 | 
-e—2(Ad)?s , (47) 


Provided the @ variation of (Ad)?, is slow compared 
to the oscillations produced by @ in the phase term, we 
may replace @ by > in the variance; integration over 0 
then yields 6(8), and 6[@—(s/g)] in (45) and (47) re- 
spectively. Integration over 8 then gives a unity re- 
flection coefficient and a ground correlation pattern 
amplitude 


e—2(Ad)’s = s/g. 


To proceed farther we must adopt a specific form for 
Mo(z). We assume a linear profile of electron density, 
beginning at height H above the ground. Then po 
= 4/1—kz, and performing the indicated integrations: 


F(6, 3) = H sec? 6 


2 
—- — (co 6 — sin? 6 — cos 6\/cos? 6 — Kz 


sin? @-cos 0 ) (42a) 
ee a 
s/cos? @ — Kz 
4 
g(0o) = 2H sec? 0) + a cos 26 (46a) 
while the variance becomes: 
k2 AN 2 a zm=cos0 /|K K222 
Ory eae f pics een 
(Ad)a , (~) Vr : (1 — Kz)3/2 
B°F2(6, 2) : cos? @— Kz 
ae ep nae gia if (44a) 
a/cos? 6 — Kz 


Evaluation of Phase Incoherence: In general this 
integral may be evaluated only in certain special cases, 
several of which we now consider. 

At vertical incidence, 0) =0, 


4 8 4H 
ro.-g-vI= Bie +z) 


4H 
[se eRe een 

=} H? -- F (42b) 

If the integration in (44a) is carried to z=1/K, a 
divergence occurs because wo =0 at this point. Since the 
first order relation (36a), for Au does not hold in this 
case, but instead u~+/AN/N, we may proceed by inte- 
grating to . 
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and then adding the contribution of the remaining seg- 
ment, 


? 


ae 


to the level of reflection. Since it is apparent that the 
major contribution comes from the region Kz~1, we 
may set =1/K in (42b). 


eae By 2 \8 
fe eS k2/AN\2 a —— jae 


N 
2=1/K (1-|AN/N 1) dz 
ab ee 
Oa ore 7 eee 
z=1/K (1-|AN/N |) AN Ay ae 
N 


=2ve el 
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Fig. 3—Ground correlation pattern of ionospherically reflected signal. 


_ Since g(0) =2[H-+(2/K) ], the scale factor in the ground 
correlation pattern is , as we might expect. For a non- 
linear profile, the electron density gradient at the level 
of reflection should be employed to determine K, since 
the major portion of the phase incoherence is intro- 
duced at this point. Fig. 3 shows a typical pattern de- 
duced from experimental measurements on ionospheric 
winds by utilizing a space-time conversion under condi- 
tions of small structure change. Since it is the signal 
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envelope rather than the rf which has been correlated 
here, a modification of the above theory is required 
which in the worst case, a Rayleigh distribution, would 
require squaring the theoretical pattern. Gaussian and 
exponential types of correlation patterns are shown for 
comparison. The best fit corresponds to parametric 
values of J~250 meters and |AN/N|~10~*. The half 
power points of the angular spectrum occur at B= + 4 

For highly oblique incidence, 00>>0,1—Kz2sin 20.~1, 
the leading terms in 


4  sin* 6) cos? 


F?(0,-2) = + H? sect 6 
K? cos? 4) — Kz 
so that 
Bly/e [AN\? eer K ae2de 
OPA MOTT acy ( = if y/cos? 0) — Kz 


2 


4 
{ a ke sin! 09 cos? 00+ H? sec* 80(cos? 80— 
1—e 


oy (42c) 


80 Rl sl 
= cos® 69 
15K N 


| { l/r (ES sec “) 
“| 4— <— —— erf | ———— 
64 GBH sec 4 ih 
H 
a5 15 —B?—— sec? Oo 
= (S + = cose 0) ¢ Y \ 


(44c) 


It is apparent that far less phase incoherence is present 
at oblique incidence than at vertical incidence. 

Finally, we consider the casé of turbulence present 
in a lower layer, in which the mean refractive index does 
not differ very much from unity, i.e, we take (AN/N)? 
as different from zero over a thickness (Az), at wo. Then: 


. Re (Az) /AN\? 
(Ad)*s,6.,a2 = = (=) Ch lee 
272 
p sec? = 


2 COS Oo 
_B 
{1 ae 
The correlation pattern in this case consists of a product 
of two exponentials, corresponding to the phase inco- 
herence introduced in the upgoing and downcoming 
waves. For the turbulent layer at height H’, and a total 
virtual height of reflection D, the scale factors have the 
values (2D —H’)/2D and H’/2D respectively. 

In all cases the frequency power spectrum may be 
obtained from the results of section 4 with “a” set equal 
to (Ad) *s-.., provided the assumption is made that the 
space-time correlation function of (AN/N) can be 
written as a product of functions of space and time in- 
dividually. This requirement arises from the integra- 
tion along a trajectory necessary to evaluate “a” in 
volume-type problems. 


(44d) 
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On the Theory of Corrugated Plane Surfaces’ 
RS. ELUIODG 


Summary—An analysis is given of an electromagnetic system 
composed of a rectangular waveguide in tandem with a corrugated 
waveguide which feeds a flat, corrugated surface of arbitrary length 
terminated by a ground plane, whose length is also arbitrary. An 
improved procedure of field determination is used which combines 
Floquet’s theorem and the variational principle, thus revealing an 
additional requirement on the corrugation geometry. Factors in- 
fluencing a match at the feed mouth, and satisfactory launching of 
the surface wave are discussed. The degree of suppression of the 
feed radiation is given in db as a function of the geometry of the 
system. Approximate radiation patterns are derived for two cases, 
(a) when the system is terminated by an infinite ground plane, and 
(b) when the system is terminated by a finite ground plane. For the 
latter case, an upper bound on the tilt angle of the main beam and a 
lower bound on its beamwidth result from an approximate theory. 
For both cases, the Hansen-Woodyard endfire relation is found to 
provide beam sharpening even when the feed radiation is considered. 
The presence of higher order surface modes, their effect, and their 
elimination are discussed. Comparison of the theory with experiment 
is reasonably good. 


INTRODUCTION 


HE THEORY OF electromagnetic waves sup- 
“ieee by corrugated conductors has been treated 

by many authors. In 1941 Slater! derived an ap- 
proximate theory for wave propagation between paral- 
lel conducting plates of infinite extent, the interior sur- 
face of one of the plates being corrugated in the trans- 
verse dimension. Goldstein?* extended this analysis to 
rectangular waveguide and Walkinshaw? applied the 
technique to circular waveguide. Cutler’ appears to 
have been the first to demonstrate that surface waves 
can exist exterior to a single corrugated conductor. By 
assuming that only the dominant surface mode was 
present, and by matching its average surface impedance 
to that of the TEM modes assumed to exist between the 
corrugations, he succeeded in deriving approximate ex- 
pressions for the phase velocity on flat and circular cor- 
rugated surfaces. By making a better approximation for 
the field distribution across the gaps between corruga- 
tions, Rotman’ obtained a slight modification of Cut- 
ler’s formulas. 


* Revised manuscript received by PGAP, January 21, 1954. Work 
reported here was performed at Hughes Aircraft Company, sponsored 
by Air Force Cambridge Research Laboratory, Cambridge, Mass., 
under Contract AF 19(604)-262, and was described in Hughes Air- 
craft Company Technical Memorandum No. 317. ; : 

+ Hughes Research and Development Labs., Culver City, Calif. 

1 J.C. Slater, “Theory of the Magnetron Oscillator,” MIT Radia- 
tion Lab. Report V-5S, pp. 1-32; August, 1941. ; a 

2H. Goldstein, “Cavity Resonators and Waveguides Containing 
Periodic Elements,” Ph.D. Thesis, MIT; 1943. . : 

3H. Goldstein, “The Theory of Corrugated Transmission Lines 
and Waveguides,” MIT Radiation Lab. Report 494, pp. 1-17; April, 
194 


4, 

4W. Walkinshaw, “Theory of Circular Corrugated Waveguide for 
Linear Accelerator,” British TRE Report T2037; August, 1946. 

5 W. Walkinshaw, “Theoretical Design of Linear Accelerator for 
Electrons,” British Proc. Phys. Soc., vol. 345; September, 1948. 

6 C. C. Cutler, “Electromagnetic Waves Guided by Corrugated 
Conducting Surfaces,” Bell Telephone Labs. Report MM-44-160-218,; 
October, 1944. oe ae 
- 7W. Rotman, “A study of single surface corrugated guides, 
Proc. I.R.E., vol. 39, pp. 952-959; August, 1951. 


All of these analyses were based on Floquet’s theo- 
rem.® In 1950, Goubau® reviewed the work of Sommer- 
feld!® pertaining to axial surface waves on cylinders and 
extended it to conductors with dielectric coats. Assum- 
ing a quasi-stationary field, he also showed that a con- 
ductor with a modified surface (e.g., a threaded conduc- 
tor) could support an axial surface mode. Lucke"™ used 
a novel approach to determine the propagation constant 
over a single, flat corrugated surface. By considering a 
finite length of surface terminated by vertical conduct- 
ing walls (which thus formed a resonator), he was able 
to get two expressions for the propagation character- 
istic—one in terms of the unknown electric field and 
the other in terms of the unknown magnetic field. The 
advantage of this method is that both formulations 
must be satisfied by the true field, so that the degree of 
approximation of any trial function may be gauged. The 
resonator method unfortunately requires the considera- 
tion of many corrugations, which makes the analysis 
cumbersome. By applying Floquet’s theorem to Lucke’s 
method, the analysis can be confined to a single cor- 
rugation, meanwhile retaining the desirable feature of 
alternative formulations. 

In the next section this modified procedure is applied 
to a corrugated rectangular waveguide to obtain the 
propagation characteristics of the allowable mode con- 
figurations. The analysis is extended to a corrugated 
parallel plate transmission line (by letting the side walls 
recede to infinity) and to a single, flat corrugated sur- 
face (by then letting the top wall also recede to infinity). 
This affords an insight into the matching characteristics 
of a corrugated waveguide feeding a corrugated surface. 

Radiation patterns of corrugated surfaces have been 
measured by the Stanford group! and by Ehrlich and 
Newkirk. Discrepancies between theoretical and ex- 
perimental patterns led Ehrlich to an experimental 
demonstration that feed radiation was contaminating 
the patterns and he was successful in devising tech- 
niques for reducing the effect of the feed. Theoretical 
support for this viewpoint will be given, together with 
curves of feed suppression as a function of the system 
geometry. It is demonstrated that by suitable design, 
the feed radiation can be reduced to reasonable propor- 


8 For a discussion of Floquet’s theorem the reader is referred to 
J. C. Slater, “Microwave Electronics,” pp. 169-177, D. Van Nostrand 
Co., Inc., New York, N. Y.; 1950. 

9G. Goubau, “Surface waves and their application to transmis- 
sion lines,” Jour. Appl. Phys., vol. 21, pp. 1119-28; November, 1950. 

10 A discussion of Sommerfeld’s analysis may be found in J. A. 
Stratton, “Electromagnetic Theory,” pp. 524-537, McGraw-Hill 
Book Co., Inc., New York, N. Y.; 1941. 

1 Second Quarterly Progress Report, “Ridge and Corrugated 
Antenna Studies,” Stanford Research Inst.; January, 1950. 

” Quarterly Progress Reports 2 through 6, “Ridge and Corru- 
gated Antenna Studies,” Stanford Research Inst.; October 1949 to 
January 1951. 1 

1M. J. Ehrlich and L. Newkirk, “Corrugated Surface An- 
tennas,” Convention Record of the I.R.E., Part 2—Antennas and 
Communications, pp. 18-33; 1953. 
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tions by some sacrifice of the match. An optimum com- 
promise can be determined by experiment. The remain- 
ing feed radiation can be phased to improve the over-all 
radiation pattern. Assuming a suitable design, the radi- 
ation pattern is then contributed to chiefly by the cor- 
rugated surface and its ground plane and is character- 
ized by an endfire main beam tilted up somewhat from 
the plane of the surface. An approximate theory is de- 
rived which permits an estimation of the pattern in 
terms of the corrugation geometry and the lengths of 
the corrugated surface and its ground plane. Curves of 
maximum tilt angle and minimum beamwidth as func- 
tions of these parameters are presented. The effect on 
the radiation pattern of the presence of higher order sur- 
face modes is found to be slight. Their presence would 
appear to be objectionable only when the corrugated 
surface is being used as a transmission line. They may 
always be suppressed simply by narrowing the teeth. 


Fig. 1—Corrugated rectangular waveguide. 


THE CORRUGATED WAVEGUIDE 


As a starting point for the analysis, assume an in- 
finitely long rectangular waveguide whose bottom wall 
is uniformly corrugated as shown in Fig. 1. It is desired 
to find the field configurations which can exist inside 
this structure. If T is the width of a tooth and G the 
width of a gap, the structure is periodic in units of 
(T+G). Hence if F(x,, y, z)e%* is the distribution of a 
field component in the plane x=x, (for OSySb, 
—h/2 S283), then F(x.+[T+G], y, z)e#t=F(x,, y, 2) 
e—#o(T+@ eit is the distribution one period further down 
the guide. 8, is a complex constant whose value depends 
on the geometry. Since the only function which satis- 
fies this requirement for all x, is e-#o* we have 
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F(a; 9; FS fee 2) ee™ (1) 


Further, f(x, y, 2) must be periodic in x in the interval 


(T+G) or 


ive] 


Ff (acy pee = gD) 2) DY ae ee 


n=—o 


(2) 


Thus each field component may be written in the form 


Da Lp naiewi sees (3) 
in which 
B nied 21n (4) 
wee 0 T 4 G 


is a complex constant to be determined by the geometry. 
This is Floquet’s theorem. 

We shall assume that the structure is to be excited 
by a Eo: mode incident from a tandem section of regu- 
lar guide. The allowable modes in the region above the 
corrugations are then hybrid, characterized by the ab- 
sence of an E, component, and given by 


ie TZ 
E, = >> An sinh an(y — b)-cos —- etn) 
n=—0o a 
a 1B, An TS . 
BE, = > Peet Gah an(y — b)* cos —- ete bn2) 
n=—0o An a 
Le ae ox ae 
H,= >, ——A, cosh a,(y — 6)-sin —- ei tfne) 
n=—0 JOUA,D a 
Tv 
tanh fe 
= a TZ 
H, = >, - sinh a,(y — b)-sin —- e7(e t Bn) 
n=—0 JOM t a 
Si fSak? : 1S 
H,= >, ——A,coshan(y — b)-cos—-ei@tn2) (5) 
n=—0 JWUAy a 
where 
aE ES Se w\2 
a, = V Bn? TKS, K= 4/® ¥. (=). (6) 
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With the origin of co-ordinates chosen as shown in 
Fig. 1, the field distribution in the gap beneath the or- 
igin may be written 


; ; eo 2s 
E, = — jou: Bo sin K(y + h) cos —-ei*t 
a 
ioe] 
marx TZ 
+ >) Bm cos ‘COS —-ef0 tr ymy 
m=1 G a 
[ova] ‘ 
mr _ Max TZ 
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mr x z 
“5 a Se Dien COS “sin ——-¢@ tym aie 
m=1 Jw a a oo 
a E = — >> A, sinh anb-e7(¢t-Bn2) (11a) 
H, = — KB, cos K(y + h)-cos—- et ome 
a co “- 
©. K2 ahs pat H = — 9) ——A, cosh and: ei(e!Bn2) (11b) 
— >>> Bm COS COS —— sete tring (7) San ached ands 
m=1 J©KLYm G a es MT Xx 
eaehiek E = — jouB, sin Kh- eit + =r cos Fae (11c) 
P mr\ af 
Vv m= ie = (8) 35 = — KB, cosKh-et#t — ss —— 
m=1 JOMY m 
Implicit in (7) is the concept that each slot is a short- Ue 
circuited waveguide of length #. The m series represents = Bm COS Gis (11d) 
that combination of TM and TE modes, all beyond 
cutoff, which together satisfy the requirement that —and 
E,=0. It is assumed that the slot width G is so small 
E= &, H = 3. (12) 


compared to the free space wavelength that all of these 
modes are attenuated to a negligible amplitude before 
they reach the bottom of the slot and hence have no 
reflected component. The one propagating mode is rep- 
resented by the standing wave terms, written separately 
in (7). Thus K and Vm, given by (8), are positive real 
quantities determined solely by the geometry. 

Since the tangential fields must be continuous across 
the boundary between the two regions, that is, across 
the plane y=O, one may equate the integrals 


+a/2 7 G4T_, 
f f Fe eg Oe 
—a/2 0 
+a/2 
= =| fé x ees. () 


in which Ei; Th, are the unknown fields expressed in 


The Fourier coefficients A, and B,, may be expressed 
in terms of integrals of the unknown electric field, i.e., 


1 Gur 
A,* = — : i Diet OR Boe) ipa 1.3)) 
(G + T) sinh a,b 
1 G 
Bo* = ——— f Exeietdy (14) 
jouG sin Khd 9 
CS Oe Eg 
Bp* = ={ Cr cos era (13) 
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and when these relations are substituted in (11b) and 
(11d) the magnetic field is given in terms of the electric 
field. Upon inserting the resulting expressions for the 
magnetic field in (10), making use of (12), and rearrang- 
ing terms, one obtains 


ee est G+ G+T PT Gian iis steal p max 
Pp ———-—-coth and f Retmede [ E*e~i®ntdy+ Xe oa Ecos ax f E* cos dx 
naw GtT ay 0 0 m=1 Ym G 0 16 
cot Kh= ; rs (16) 
f Edx al E*dx 
0 0 
By asimilar procedure, (12) can be expressed entirely in 
terms of the unknown field H, giving 
G+T G+ 
> pan mee oc ant aii ed Heinzdx: H* cthedt Dy - —- ef, Heme logon de ax f fairy ae 
1 _ n= G+T wre 0 
cot Kh 


terms of (5) and é,, 5, are the unknown fields expressed 
in terms of (7). All terms of the integrands of (9) con- 
tain z only in the factor cos 7z/a. Hence (9) reduces to 


G+T G 
f EH*dx = f Esxc*dx (10) 
0 0 


ery H*dx 


Equations (16) and (17) are suggestive of Schwinger’s 
variational form but investigation discloses that only 
(17) is stationary about the true fields.* However, both 
must be satisfied by the true fields and from this fact 

“4 R, S. Elliott, “On the Theory of Corrugated Plane Surfaces,” 


Hughes Research and Development Labs., Technical Memorandum 
No. 317; October, 1953. (Revised.) 
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one can derive useful information. To see how good an 
approximation the fundamental mode is, we insert the 


trial function E=—Apo sinh a,b-e%*%” in (16) and 
obtain the asymptotic formula 
Gre es 
cot Kh = -— coth ab (18) 
Qo 


as G+T-—0. Similarly, when the companion trial func- 
tion 


7 2 


H=-— A, cosh ayb- e7(@ For) 


JOU 


is inserted in (17), the asymptotic formula 


K 
-—-coth a,b 
GT a 


cot Kh = 


(19) 


results when G+T7-—0. Previous work by the Stanford 
group on a similar problem" indicates that these form- 
ulas are good representations for (A,/G+T) 210. 

It is interesting to note that only for 7=O are the 
two formulas the same. This is reasonable when one 
recalls that Ei,an=O over each tooth and the trial func- 
tion used to obtain (18) does not satisfy this require- 
ment if T>O. Hence, we conclude that only if (a) the 
number of corrugations/wavelength is large, and (b) 
tooth width/gap width is small, does the field distribu- 
tion above the corrugations consist essentially of the 
fundamental mode. Point (b) seems previously to have 
been overlooked and it can have important bearing on 
the impedance concept when corrugated surfaces are 
used as transmission lines. This shall be discussed fur- 
ther in a later section. 

The question still remains as to which formula, (18) 
or (19), is more accurate when T>O. One would sus- 
pect (19) is, because of its stationary character, and be- 
cause of the severe requirement on Eitan. Support for 
this belief arises when any more general trial function 
is substituted in (16) and (17). The variation in (16) is 
always greater than the corresponding variation in (17). 

Henceforth we shall assume that there are at least ten 
corrugations per wavelength and that the teeth are nar- 
row compared to the gaps. Then to a good approxima- 
tion the field above the corrugations is given by 
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(gy Bo? so KA 
K 
cot Kh = -—-coth a,b 
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2 
Kee 4/ # = (=). (20) 
a 

We notice that if OX KhS(a/2) then 8, and a, are both 
real and 8. >K so that wave propagation is slower than 
it would be in uncorrugated guide. 

If the side walls are allowed to recede to infinity, the 
solution for a corrugated parallel plate transmission line 
is obtained. Namely, 


Ex = Agisinh ao(y i= bette! Po”) 


iE = eee cosh aly = b)e7(e # Box) 
Qo 
Jue 
H, = — A, cosh a(y — bet For) 
Qo 
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GET Bol) St 
If the spacing 6 becomes indefinitely large, the solu- 
tion for a single flat corrugated surface results: 


Nabe = Ae t— 7jBot—agy 
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in which 


(22) 


Thus we conclude that if the side and top walls of a cor- 
rugated waveguide are gradually flared out and then 
terminated, a surface wave may be satisfactorily 
launched on the extended bottom wall with the expecta- 
tion of a good match. If at the other end of the cor- 
rugated waveguide the depth of corrugations is gradu- 
ally tapered to zero, a good match to regular guide can 
be achieved” and the resulting system can be efficiently 
excited by a TE, mode.% 


‘ 


RADIATION FROM CORRUGATED SURFACES 


A problem of considerable practical interest is the 
computation of the radiation pattern of the system of 


® As with any horn, there are space limitations to such a system 
and for wide surfaces, line feeds (such as pill boxes, hog horns, etc.) 


are more desirable. In such cases, flaring of th is sti 
ee eee es, flaring of the top wall is still helpful 
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GROUND 
PLANE 


Fig. 2—Corrugated surface antenna. 


Fig. 2. We shall assume all the requirements for a good 
match are met, i.e., tapered corrugations in the wave 
guide and gentle flare angles for the horn to a large 
aperture bxw. It is then reasonable to expect that most 
of the power incident in a TE»; mode is transformed to 
a surface mode as given by (22). 

Since the chief effect of the finite width w is to alter 
the horizontal beamwidth, we shall infer the solution of 
this problem from the similar but simpler problem of an 
infinitely wide parallel plate transmission line whose 
lower plate is corrugated and extended out to form a 
single surface, as shown in Fig. 3. For the present, we 
shall assume the surface to be terminated by an infinite 
ground plane, deferring to the next section a discussion 
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implies that the currents which leak back over the out- 
side of the upper plate are negligible. If f<l it may be 
ignored. If not, the effective length J’ can be taken as 
some reasonable compromise, such as /’ =1+#f.) 

Lucke has shown" that for f=O the reflection coeffi- 
cient for a surface wave of the form (22), incident at the 
junction of the corrugated surface and its infinite ground 
plane, is given in magnitude by 


Bo =. k 
Bo 


[R|= (23) 


We shall see shortly that maximum gain results for the 
corrugated surface of Fig. 3 if (8,—k)l=7 (the Hansen- 
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Fig. 3—Simplified antenna system. 


of the effect of a finite ground plane. The radiation pat- 
tern of the system of Fig. 3 is contributed to by the 
secondary source distribution across the mouth of the 
aperture, the current distribution in the corrugated 
section of length /, and the current distribution in the 
infinite ground plane. (The assumption of a good match 


Woodyard relation). Hence for practical systems B,/k — 1 
is small and likewise | R| is small. The tapered section f 
tends to reduce | R| still further so we shall assume no 
reflection at all. 


16 Third Quarterly Progress Report, “Ridge and Corrugated An- 
tenna Studies,” Stanford Research Inst.; April, 1950. 
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The field distribution over the corrugated surface is 
then given by (22). If this field is terminated by the 
electric and magnetic current sheets 


= = JWE ; 
j-= — 4, Ages) (24) 
A,* 
= = : 
Jp Ug Aer (25) 


the corrugated surface may be removed and the ground 
plane extended back to the mouth of the feed. The ef- 
fect of the ground plane can then be accounted for by 
the method of images. 

Referring to (21), the field at the mouth of the feed 
can be terminated by the sheets 


= — Jue 
J-=- Su, : - A, cosh ak (y — b)ei(@t-Bo" 20) (26) 
Ao 
> = ipF = 
Im = — Uz Se A, cosh aF(y — b)ei@t Bo 0) (27) 
en 


in which x, is the «-co-ordinate of the aperture. (The 
superscripts s and F are employed as mnemonic devices 
for the surface and feed respectively.) The images of 
(24)—(27) are 
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Qo 
so that the original system is equivalent to a magnetic 


current sheet 
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occupying the same position as the corrugated surface, 
and a double sheet 
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extending from —b to +0 in the plane of the feed 
mouth («=x,). 
The radiation patterns of these sources are 
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_ gilwt—kpotr/4) gi(n LIN) (Bo. | k—c0s0) 
k sin 6-sin [kb sin 6| + a" sinh a” 
(asf)? + k?.sin? 6 


(32) 


in which the center of the corrugated surface has been 
chosen as origin (making «,=—I/2). The corrugated 
surface is seen to give a conventional endfire pattern 
which can be maximized by setting 


(B.% — k)l = 7m. 


This is the Hansen-Woodyard relation mentioned pre- 
viously. 

Although there is considerable turbulence in the re- 
gion of the mouth, if the guide and surface modes are 
extrapolated to the position x = —//2, we can write 


(33) 
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in which C and wy are positive real numbers which de- 
pend on the corrugation geometry and the mouth 
height, 6, but not on length of corrugated surface, /, 
since it is assumed there are no reflections. The quantity 


C? | ae 
sinh? eb 


can be found by equating the power in the guide and 
surface modes, which gives 


(35) 


(36) 


2 Bo° (ask \3 1 
dal GE (=) a¥’b + sinh a,¥b-cosh obo 
Both fields (31) and (32) have their maxima in the 


direction 6=0 degrees. Making the substitutions (33) 
and (34), the ratio of these maxima is 


\ ‘HA 2B =k) 
--——_— (37) 
ie aCe 
and the corresponding power ratio is 
Ps 4 4]?(B.F — k)? 
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This relationship is plotted in Fig. 4 for several lengths 
of corrugated surface. It is seen that the mouth must be 
made quite small to effect satisfactory feed suppression. 
This is going contrary to the requirements for a good 
match, and the analysis breaks down when the match is 
poor enough that it is no longer proper to equate the 
powers in the two modes. Hence, some optimum height 
b exists, which for practical systems probably will not 
yield more than 10 to 15 db of suppression. This opti- 
mum height will be a function of 8,* and J and can best 
be determined by experiment. 

It is observed from (37) that ify is small (which seems 
reasonable) the two fields are essentially out of phase in 
the forward direction. Since the feed pattern is broad, 
this tends to sharpen the main beam. As an illustration 
of this, Fig. 5 shows the experimental pattern of a 7.33X 
corrugated surface as measured by Ehrlich. 8,* was ad- 
justed to give maximum gain, ie., (8.%—k)l=7. For 


NOTE: 


THEORETICAL UPPER BOUND 
ON THE SUPPRESSION OF 
FEED RADIATION FOR THE 
SYSTEM OF FIGURE 3. 


Fig. 4—Feed suppression in db. 


comparison, the theoretical patterns for the corrugated 
surface alone, and for the surface plus feed with y=O, 
m are plotted. It is seen that the case y =O corresponds 
most closely to experiment. Thus it is wise to use the 
Hansen-Woodyard relation even when the feed radia- 
tion is considered, for the phasing is then proper to 
provide additional beam sharpening. This will also be 


seen to be true in the presence of finite ground planes. © 


THE EFFECT OF A FINITE GROUND PLANE 


Thus far the corrugated surface has been assumed to 
be terminated by an infinite ground plane which per- 
mits the use of the image principle and greatly simpli- 
fies the analysis. A reasonable approximation to the 
pattern for the case of a finite ground plane can be ob- 
tained by assuming the dominant mode to be present 
over the corrugations and by assuming the same current 
distribution to exist in the finite ground plane as would 
exist in the same portion of an infinite ground plane. 

For the case of no ground plane at all, this implies 
the radiation pattern arises essentially from the double 
sheet (24) and (25). (We assume feed radiation is suffi- 
ciently suppressed as to be only a minor perturbation on 
results following.) This radiation pattern is given by 


welA / Z 
4a,? Tv Rpo 


H, = — (a* + jk sin 6) ei(wt—kpot m/4) 
_ wT B® 
ee ; — cos @ 
ee 19" cea (39) 


and is an approximation which is valid only for (8,*—k)/ 
small. A plot of (39) is characterized by a main beam 
tilted up from the endfire position and slightly broad- 
ened with respect to the infinite ground plane case. 
These results are consistent with experiment. The tilt 


angle of the main beam is given by 
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Fig. 5—Theoretical and experimental radiation patterns for a 
corrugated surface of length 7.33 excited by a waveguide. 
This tilt angle is a decreasing function of 6,°/k for all 
lengths / and thus has an upper bound for 8,*/k=1. 
This is fortunate, for the approximation is most accu- 
rate for this minimum value of 8,°/k. A plot of the upper 

bound is shown in Fig. 6, on following page. 

To see what happens as ground plane is added, we 
first must find the current distribution in the infinite 
ground plane. Using the source system (28), we find 


weA Oye: Z 
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in which u=k(x—<x’). When the Hansen-Woodyard 
relation, (8.°—k)l=m7, is satisfied, the phase change of 
the integrand of (42) over the interval k(x—1/2)Su 
<k(x+1/2) is approximately 7 radians for all «. Thus 
the integrated error accruing from the substitution 
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Equation (45) may be solved in terms of Fresnel inte- 
grals through the substitution 
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Fig. 7—Theoretical upper bound on the main beam-tilt angle of the 
radiation pattern for a flat corrugated surface as a function of 
the lengths of the surface and its ground plane. 
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Since the value of J; is also known at x =1/2, the current 
distribution in the intervening half-wavelength can be 
inferred by extrapolation. Equation (47) indicates that 
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Fig. 8—Comparison between theory and experiment for two corrgu- 
gated surfaces obeying the Hansen-Woodyard relation, one 7.33 
long and the other 4) long. 


‘ 


for a given x, the current density is a decreasing function 
of 6.°/k. Thus for a finite ground plane of length d, with 
the assumed current distribution (47), the larger the 
value of 8,°/k, the smaller the end disturbance and the 
smaller the tilt angle. Therefore the tilt angle is a de- 
creasing function of d for all 8,°/k, with the rate of de- 
crease least for 6.°/k =1. But for the case 6.?/k=1, the 
corrugated surface of length J, and its ground plane of 
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Fig. 9—Geometrical optics formulation. 


length d are electrically indistinguishable, and Fig. 6 
becomes the upper bound when (J+d)/\ is substituted 
for //X as the abscissa scale. This information can be 
plotted in an alternative way as shown in Fig. 7, left. 

To check this theory, a series of ground planes of 
length 0A to 19\ were added to a 4 corrugated surface 
and the tilt angles measured. In a second experiment, 
ground planes ranging from 2) to 8A were attached toa 
7.33X corrugated surface and the measurements re- 
peated. The results are shown in Fig. 8, page 78. The 
agreement with theory is fair, with the greatest de- 
parture occurring for long ground planes. 

A satisfactory picture of the general nature of the 
four-quadrant radiation pattern can be gleaned from a 
geometrical optics argument. 

Referring to Fig. 9, (shown above) if we neglect the 
feed radiation and assume that the current distribution 
in the ground plane does not radiate in the backward 
direction, the radiation field along the half-plane A —A’ 
is given by (39) and if this field is terminated by the 
proper double sheet the radiation field in quadrants II, 
III, and IV is approximately determined by this double 
sheet in the presence of the ground plane as an obstacle. 
The amplitude distribution of the source system along 
A-—A’ is suggested by the curve f(s) and serves to ex- 
plain many of the features of the pattern. 

For an angle 6; in quadrant II, the radiation is com- 
puted from the sources along A—A’ plus the image 
sources from B to A. As 6; approaches 90 degrees, the 
radiation approaches the value computed for an in- 
finite ground plane, and throughout quadrant I we as- 


sume that the radiation pattern does correspond to this 
infinite case. As 6, approaches 0 degrees, fewer of the 
image sources contribute, and in the position 62, the 
field strength is down to one-half the value found for 
the infinite case. At any angle 6,, the radiation ap- 
proaches more nearly to the infinite case as the ground 
plane is lengthened. 

For an angle 63 in quadrant III, only the sources from 
C to A’ contribute and as 6; approaches —90 degrees, 
field decreases to zero, oscillating slightly as changing 
source system phases in and out. Under this geometrical 
argument, no energy is found in quadrant IV. 

Using the above model, a sketchy picture of the radia- 
tion pattern may be synthesized as follows: For the 
lengths / and d of surface and ground plane being con- 
sidered, determine from Figs. 6 and 7 the approximate 
tilt angle, 67, of the main beam. Use the field distribu- 
tion (31) for the region 6r $6180 degrees, assuming 
the null positions are undisturbed. For the region 
—90 degrees S067 sketch in a smoothly decreasing . 
field approximately 6 db down at 6=0 degrees. Result 
for a 7.33 surface with a $A ground plane is in Fig. 10, 
page 80. This pattern agrees in its general shape with 
measured patterns, two examples of which are given in 
Fig. 11, on page 81. 


BEAMWIDTH 


The geometrical optics picture just outlined suggests 
that for a corrugated surface of given length “/,” as the 
ground plane length “d” is increased, the main beam not 
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Fig. 10—Empirical radiation pattern for a 7.33\ corrugated surface with a /2 ground plane. 


only tilts down closer to the plane, but the beamwidth 
narrows, approaching the infinite plane value as a limit. 
Plotted in Fig. 12, on following page, are the results of 
beamwidth measurements for a 7.33 surface and a va- 
riety of ground planes. As can be observed, the pattern 
actually narrows as the endfire position is approached. 

A lower bound on the beamwidth can thus be found 
from (31). A plot of this lower bound as a function of 
1/X is shown in Fig. 13 for (6,—k)]=7 on the following 
page. 

THE EFFECT OF FINITE TootH WIDTH 


We have already observed that unless the teeth are 
extremely narrow, the dominant mode cannot ade- 
quately satisfy (16) and (17) and that the wider the 
teeth, the more the higher order modes are excited. Re- 
ferring to (4), if there are twenty corrugations per wave- 
length, 

By 
— ~~ 


leet SANT aA etc. 
k k 


and if these modes are excited, they will have radiation 
patterns expressed in the form of (39), with appropriate 
substitutions for a, B, and A. The factor 


(48) 


(an® + jk sin 6) (49) 
is practically insensitive to 6. The factor 
Lae | [ope 
Se = iCOS | 
Ge: 
(50) 


gives essentially the same pattern for all higher modes, 
i.e., 21/A equal lobes with the nulls occurring at different 
angles for different modes. The fields of these higher 
modes will have phase relationships which depend on 
the complex constants A». 

An effort was made to excite these modes by con- 
structing a corrugated surface with twenty teeth per 
wavelength and teeth three times as wide as the gaps. 
Its radiation pattern was compared with that of a sim- 
ilar surface in which the gaps were three times as wide 
as the teeth. The slot depths of the two surfaces differed 
by the amount necessary to establish the same value of 
B.°/k, a value chosen to satisfy the Hansen-Woodyard 
relation. Two typical patterns are shown for compari- 
son in Fig. 11. No essential differences were observed. 
Though not reported here, this experiment was per- 
formed over a range of frequencies, for several ground 
plane lengths, including an effectively “infinite” ground 
plane, and for tooth width/gap width ratios of 3:1, 1:1, 
and 1:3. In no case were there any significant differences 
in the patterns. The conclusion appears that the higher 
order modes are excited when the teeth are wide but 
that their amplitudes are comparable and phases random 
so that their net effect on the radiation pattern is small. 

It should be emphasized that this does not mean that 
the effect of the presence of higher order modes is 
negligible when the corrugated surface is being used as 
a transmission line. In that case, the impedance concept 
is vastly complicated by their presence and it would 
seem advisable to use narrow teeth in an effort to obtain 
an uncontaminated dominant mode. 
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Fig, 11—Experimental radiation patterns for two corrugated surfaces 7.33 long, one having a tooth to gap ratio of 3:1 
and the other a tooth to gap ratio of 1:3, each terminated by a A/2 ground plane. 
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Fig. 12—Beamwidth versus ground plane length for a 
7.33 corrugated surface. 


CONCLUSIONS 


The analysis of a flat corrugated surface appears to 
be strengthened by a Floquet-Lucke method of field de- 


_ termination. Consideration of feed radiation gives bet- 


ter correlation between theory and experiment. Effect of 


_a finite ground plane can be explained approximately by 


classic assumption of current distribution. Effect of fi- 


~ nite tooth width, though it complicates transmission line 


analysis, apparently has no influence on radiation pat- 
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Fig. 13—Theoretical lower bound on beamwidth as a function of 
length of corrugated surface when Hansen-Woodyard relation is 
obeyed. 


tern when corrugated surface is used as an antenna. 
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Meteoric Radio Echoes* 
LAURENCE A, MANNING} 


Summary—Study of meteoric effects by radio methods has re- 
sulted in gain to three fields—astronomy, upper atmosphere physics, 
and radio propagation. In each of these fields progress is surveyed, 
and some of the more important unsolved problems are pointed out. 


INTRODUCTION 


ETEORS have been seen for centuries, photo- 
| \ | graphed for decades. Only since World War II 
have they been extensively studied by radio 
echo methods. In these past few years, though, the rate 
at which facts have been gathered about meteoric ef- 
fects has been the greatest in history. It is already time 
to appraise how far we have gone, and to ask what has 
been neglected in our haste. This paper, then, is a sum- 
mary of progress in radio meteoric detection; but it tries 
to be more—it tries to classify work done to date, so 
as to make evident what is yet undone. 

Meteoric studies can be placed in one of three groups. 
In group one are studies of the meteoric particles, 
treated as astronomical bodies; size, mass, velocity, or- 
bit, composition are unknown. In group two are studies 
of the upper atmosphere; meteors are tools used to 
tickle the air, excite a response. In group three are 
studies of the influence of meteoric reflections on the 
propagation of radio signals. As a guide to the discus- 
sions which follow, we may list under these three groups 
the main fields of past and foreseeable future interest to 
the radio worker: 

(a) Astronomical properties of meteors—1. Radiants. 
2. Velocities. 3. Size and rate of arrival. 

(b) Meteoric effects as a tool in finding atmospheric 
properties—1. Ion production, diffusion, recombina- 
tion, composition. 2. Winds. 

(c) Meteoric contributions to propagation—1. Ob- 
liquely propagated high frequency signals. 2. Ambient 
ionization and low frequency reflections. 


In each of the fields above, some work has been done. 
How much? 


Stupy oF METEORIC RADIANTS 


Meteoric radiants are found using only the simplest 
property of meteoric ionization trails—their shape. A 
typical meteor leaves a trail of ionization along its path 
that is about 25 kilometers long and but a few meters in 
diameter. Pierce, working at Harvard, was the first to 
point out that the very small ratio of diameter to length 
of these trails requires that a normal from the observa- 
tory strike the ionized part of the path, if a strong echo 
is to be received.! The reason for this behavior may be 

* Manuscript received by PGAP, September 24, 1953. Jointly 
supported by the U. S. Army Signal Corps, the U. S. Air Force, and 
the U.S. Navy (Office of Naval Research). 


} Electronics Research Lab., Dept. of Electrical Engineering, 
Stanford Univ., Stanford, Calif. 


1J. A. Pierce, “Abnormal ionization in the E-region of the 
ionosphere,” Proc. I.R.E., vol. 26, pp. 892-902; 1938. 


seen from the geometry of Fig. 1. The circular arcs of 
the figure represent equiphase surfaces of the incident 
wave. Assume that the surfaces are one-quarter of a 
wavelength apart; energy returned to the transmitter 
from point Po will clearly differ in phase by one-half 
cycle from that returned from P;, but be in phase again 
with that from Ps, Ps, Pon, etc. The part of the trail 
from P_, to P, is called the first Fresnel zone. It is 
longer and returns a stronger signal than does the part 
from P,; to Pe, or any of the higher order zones. Since 
the phase of signals returned from the higher zones 
alternates in sign, their contributions largely annul 
themselves. It is only when the first zone (centered 
about the normal reflection point) is ionized, that a 
strong signal is returned; conversely the presence of a 
reflected signal is an indication of normal reflection in 
the high majority of cases. The main exceptions occur 
for large meteors, and the lower frequencies; they seem 
to go with ionization distributed in other than the usual 
linear model. 


Fig. 1—The normal reflection condition for a meteoric ionization 
column. The circular arcs represent spherical wavefronts spaced 
by one-quarter wavelength. The scale is of course exaggerated; 
the first Fresnel zone P_;—P; should be only about 1 per cent of 
the range OPo. 


A number of ways have been developed to use the 
normal reflection property in finding radiants. As is 
indicated in Fig. 2, it is assumed that an echo will occur 
only when the reflection point is within a height region 
of relatively narrow bounds. Let hk; and hz be those 
bounds. The intersection of the plane » perpendicular 
to the radiant direction, with the region between 
heights /; and he, is drawn as the vertical hatched band 
in the plan view of the figure. All reflection points must 
lie in this band for a given radiant position, and so the 
position and orientation of this band serve to fix the co- 
ordinates of the radiant. 

Clegg has made use of the localization of the reflecting 
regions to find radiants by setting up a highly directive 
antenna beam.’ This beam is shown in the plan view of 


_ ?J. A. Clegg, “Determination of meteor radiants by observa- 
oe a mee echoes from meteor trails,” Phil. Mag., vol. 39, pp. 577- 
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Fig. 2. Over a span of time, earth rotation causes a 
change in radiant direction, and so of the position of 
the reflection point band. By measuring range to that 
part of the band in the beam versus time, Clegg is able 
to locate the meteoric radiants. He has reported preci- 
sion of about two degrees in the description of radiant 
direction. Comparable beam widths are not needed in 
the directive antennas. 
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Fig. 2—The location of reflection points for meteors of radiant R. 
Ionization is assumed formed only between heights H, and Ap. 
The reflection point must lie in the plane P perpendicular to the 
trail direction. 


It is clear that many other arrangements of antenna 
beams could be used to get the same result—the loca- 
tion of the band of possible reflection points. For exam- 
ple, an automatic scanning direction finder has been 
used in conjunction with range information at Stan- 
ford; in this way the position of the band can be plotted 
directly, without awaiting earth rotation. Variations of 
these techniques have been used by other workers.’ 

The methods of radiant determination which have 
been described have one main drawback; they are based 
upon analysis of a large number of meteoric echoes. At- 
tempts have also been made to find radiants of single 
meteors by radio methods. The writer suggested the use 
of a triangulation technique in an early discussion of the 
problem.‘ In this way, it was shown that if the ranges 

3 J. S. Hey and G. D. Stewart, “Derivation of meteor stream radi- 


dio reflection methods,” Nature, vol. 158, p. 481; 1946. 
ge PA Momadine, “The theory of the radio detection of meteors,” 
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and relative times of echo arrival were determined for 
three ground stations, the perpendicularity of the re- 
flection process could be used to locate both the direc- 
tion in space from which the particle came, and the point 
on the ground met by the extension of the trail. The 
differences in echo starting times are a measure of the 
relative distances along the trail to the intersections of 
the normal rays from the observing stations; this fol- 
lows because a strong echo is not received until ioniza- 
tion is created at the perpendicular point. Application of 
the technique has been made by Millman and McKin- 
ley, but with the use of range-time plots to determine 
when the particle and its wake of ionization reach the 
nearest (the normal) point on the trail.> The triangula- 
tion technique has the disadvantage, in addition to the 
need for three isolated installations, of a relatively low 
rate of coincident meteor detection. There are not many 
orientations of a line 25 km long and 100 km in altitude 
which permit a normal from each of three stations sepa- 
rated the order of 40 km. 

At present it seems that the most promising lines for 
future development of methods for finding the radiants 
of single meteors lie in the utilization of either fading or 
polarization effects. The use of fading properties would 
depend on the observed fact that certain trails break 
up, Causing an uneven reflection coefficient along the 
length of the column. The simple Fresnel theory does 
not then apply; wind drift causes the signals of irregular 
strength scattered from parts of the trail of different 
range to beat together to form a fluctuating received 
signal. It can be shown that under certain conditions, 
by noting the phase of the fading envelope at three sta- 
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Fig. 3—Reflected signal strength versus time for the two character- 
istic polarizations of a typical meteor echo, observed on 23.1 mc 
at Stanford University by Van Valkenberg, using the two helix 
technique. 

tions separated a kilometer or so apart, the horizontal 

projection of the trail direction can be obtained. In 

combination with horizontal and vertical elevation angle 
recording, it may be possible to use this effect for radiant 
fixes. The use of polarization properties for finding radi- 
ants also shows promise. The signals reflected when the 
transmitted electric field vectors are parallel and per- 
pendicular to the trail are not equal during the early 
phase of an echo. Fig. 3 shows the amplitudes ob- 


5 P. M. Millman and D. W. R. McKinley, “Three-station radar 
and visual triangulation of meteors,” Sky and Telescope, vol. 8, p. 
114; 1949, 


84 LR.E. TRANSACTIONS—ANTENNAS AND PROPAGATION 


(a) SUM OF ELEMENTS 


April 


FRESNEL ZONE 


(b) coRNU SPIRAL 


Fig. 4—The way in which the vector signals returned from the successive elementary lengths of the trail combine at the receiver: (a) 
Each vector represents the signal returned from a unit length of the trail; the phase of the vector depends upon the two-way path- 
length from the element to the receiver; (b) As the element length approaches zero, a Cornu spiral is obtained in the limit. Arc length 


corresponds to trail length. 


served for the two polarization components of a typical 
echo. Again in combination with directional informa- 
tion, the polarization effect may be useful for finding 
individual radiants. However phase shifts may exist be- 
tween the parallel and perpendicular reflection coeffi- 
cients; the direction of the major axis of the reflection 
coefficient ellipse would then be skewed with respect to 
the trail. This effect would need to be looked into.®7 


Stupy oF METEORIC VELOCITIES 


Three ways of measuring velocities are of importance. 
In order of their usefulness, they are—1. The Doppler 
meteor-whistle technique; 2. The diffraction technique; 
and 3. The change of range technique. All three methods 
are based upon the same properties of the ionization 
columns. In each case it is assumed that the ionization 
is created at once upon the arrival of the meteor. Thus, 
at any moment during the formation of the trail the 
ionization exists along a half-infinite line. : 

In the change of range technique, the slant range to 
the nearest part of the ionization column is plotted 
versus time during the column formation.® From analy- 
sis of the observed hyperbolic relationship velocity is 
deduced.® The main trouble is the small percentage of 
echoes for which it is possible to follow the range varia- 
tion very far from the perpendicular reflection point. 
When the trail is in formation, only the short weakly 
reflecting high order Fresnel zones are effective. To get 
enough range difference to resolve the hyperbolic con- 
stants, the echo must be observed for a long period of 
time, with the resulting high order Fresnel zones and 


SJ. A. Clegg and R. L. Closs, “Plasma oscillations in meteor 
trails,” Phil. Mag., vol. 39, p. 577; 1948. 

7M, E. Van Valkenburg, “Polarization and Fading Studies of 
Meteoric Echoes,” Ph.D. Dissertation, Stanford Univ., Stanford, 
Calif.; 1952. 

8 J. S. Hey, “The Giacobinid meteor shower, 1946,” Nature, vol. 
159, pp. 119-121; 1947. 

9D. W. R. McKinley and P. M. Millman, “Determination of the 
elements of meteor paths from radar observations,” Canad. Jour. 
Res., vol. A27, p. 53; 1949. 


weak echoes. Associated with the small rate of detection 
of usable meteors is a relatively low accuracy on all ex- 
cept the largest echoes. 

The Doppler whistle and the diffraction techniques 
are very close cousins theoretically.!° However, they 
involve different observing techniques, and lead to 
rather different sensitivities and accuracies. It will be 


_well to analyze them both from the same general view- 


point. Consider the trail to be divided into many ele- 
ments of differential length. Assume also a uniform and 
not too dense ionization distribution along the trail. 
Then each of the elemental lengths of trail will return a 
component of reflected signal. It may be described by a 
vector with magnitude proportional to the length of the 
element, and with a relative phase determined by its 
position in one of the Fresnel zones. The total signal 
strength returned from the trail at any time during its 
formation is then found by adding up the vectors cor- 
responding to each element of trail, all in correct phase 
relation. In Fig. 4(a) the manner in which these ele- 
mentary signals add up is shown. Fig. 4(b) shows the 
theoretical Cornu spiral which is obtained by letting the 
length of the elements approach zero. In Fig. 5(a) is 
shown the way in which the amplitude of the reflected 
signal vector increases with time. The oscillation of 
amplitude following the initial build-up is the “diffrac- 
tion fluctuation” utilized to measure velocity in the 
diffraction technique.!! As can be seen by comparing 
Figs. 4 and 5, the period of the fluctuations is related to 
the rate at which the trail is formed beyond the per- 
pendicular point, and corresponds to the revolutions in 
the upper turns of the Cornu spiral. 

In the Doppler meteor-whistle technique, a fixed sig- 


10 L. A. Manning, O. G. Villard, Jr., and A. M. Peterson, “Double 
Doppler study of meteoric echoes,” Jour. Geophys. Res., vol. 57, pp. 
387-403; 1952. 

"C.D. Ellyett and J. G. Davies, “Velocity of meteors measured 


by diffraction of radio waves from trails during formation,” Nature, 
vol. 161, pp. 596-597; 1948. 
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nal of large amplitude is added to the received signal, THE RATES OF METEORIC ARRIVAL 
so that a large signal is in the receiver at all times.!2 
When an echo is received with a beating voltage pres- 
ent, it is possible to identify both spiral arms of Fig. 
4(b) as an amplitude fluctuation. The first oscillation is 
used the most in the Doppler system. Among the merits 
of this method note that cross-correlation or coherent 
detection is naturally used, thereby improving the weak 
signal sensitivity. 


To the astronomer, the rate of occurrence and dis- 
tribution in size of meteoric particles is of great impor- 
tance. But it is not as easy to count meteors of a given 
size as it is to measure speed or direction. The trouble 
is twofold: observed counts are less than the true rate 
because of aspect, range and other effects; particle size 
can be found from echo amplitude or duration only with 
the aid of very shaky theory. 

In his 1938 paper,! Pierce made an important con- 
tribution to the theory of meteor counting. Based upon 
the aspect sensitivity inherent in his normal reflection 
hypothesis, he was able to show that only a fraction of 
those meteors entering the atmosphere within normally 
detectable ranges (from a random radiant) would ac- 
7 —— tually be detected. For example, trails directed straight 
Fig. 5—The theoretical variation of received signal strength during at the observer do not have a normal for reflection with- 

trail formation. In (b) a beating voltage of magnitude B has in the 80-120 km height range, and so will not usually 

ee ee medl cignalnf 49): give an echo. Manning extended Pierce’s calculations to 


Ellyett and Davies first described the diffraction tech- the case of shower meteors from point radiants.* The 
nique for velocity determination.!! Following a sugges- Same type of effect occurs; commonly the observed 
tion of Herlofson they were the first to recognize the meteor count using nondirective putennas may be only 
importance of the Fresnel approach. The observation of 10 per cent of the true rate of arrival. It is therefore 
the Doppler effect had been made many years before "ecessary to correct the observed rate. The actual size 
by Chamanlal and Venkatamaran." It was later inde- of the correction is hard to fix, since it depends on such 
pendently noted by Villard.'4 The first theoretical treat- usually unknown quantities as radiant distribution and 
ments of the Doppler whistles were by Appleton and trail length. Manning, Villard and Peterson have 
Naismith*® and by Manning.‘ Later Manning, Villard, evolved a method for getting values of trail length; if 
and Peterson verified the accuracy of the technique ex- applied more generally it would reduce some of the un- 
perimentally, using the Perseid shower as a source of certainty.!” If the rate is to be found for shower meteors 
meteors of beketc velocity.!2 McKinley has measured from a known radiant, directive antennas may be used. 
the velocities of thousands of meteors by the Doppler Range data will then reduce much of the uncertainty 
method. He firmly established the usefulness of the about the fraction detected in the most favored direc- 

ea i H : tion. 
technique by showing that no appreciable fraction of , wee 
radio meteors have velocities corresponding to hyper- aa yet reich erate an amie aiden 

i its.16 therefore members of our own With a given echo. itis ntly that 
ies : ae pee theory of radio reflections from meteoric ionization trails 
so F . ° . . 

In Ho to higher sensitivity, the Doppler tech- has been available, from which even the relationship 
nique has an advantage accruing Mea the experimental between ionization production and echo strength could 

al that the diffraction fluctuation of many ob- be reliably inferred.!* No reliable theory 1s yet available 
Be meteor echoes is weaker, and contains fewer 1° relate meteoric mass and ionization production. When 
cycles, than the preceding ede whistle. As many as such a theory is at hand, however, it seems certain that 
three ed four hundred cycles of down-whistle oscillation radio methods will be able to augment with precision 
have been recorded; the up-whistle or “diffraction ef- the existing knowledge, particularly with respect to the 

’ . . . . 

fect” rarely can be followed so long. The explanation of distributions of the smaller meteoric masses. 
this effect has never been properly sought, although it is Paocrasrs IN Tae LONIZATION (RAID 
known that the ratio of down- to up-whistle amplitude 
increases with echo size. 


DIFFRACTION FLUCTUATION 


Meteors are of use in finding the physical properties 
of the upper atmosphere because they produce a tran- 
#1, A. Manning, O. G. Villard, Jr., and A. M. Peterson, “Radio sient increase in ionization; the resulting time response 


i igati ic hei ities,” Jour. Appl. j d 
Fine val 20, te ree ay all fl Ras er HPP of the medium may be easily observed. The details of 


18 C, Chamanlal and K. Venkatamaran, B a4 ai cig aie the ionizing process are not understood; however, the 
Dover ea sec al Picbictaa game way in which the ionization is distributed about the 
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meteoric path, and the type of reflection which then 
results, are subjects that are finally yielding to experi- 
ment and analysis. The early ideas about the amount 
of ionization were those of Skellet!® and of Pierce.”° 
They postulated that the ionization surrounding the 
meteoric path extended, at least in the case of large 
meteors, to radial distances of up to a kilometer. These 
ideas of size were based in part on analogy with the 
glowing trains which are sometimes seen visually and 
telescopically.?! Justification for the amount of ioniza- 
tion was given by assuming that nearly all of the energy 
of the meteor went to ionize the surrounding air. Most 
experiments at high frequencies are not consistent with 
the assumption of columns of the breadth assumed in 
the early work. Perhaps at the frequencies of 3.5 mc or 
so used by Pierce it may be possible to detect low 
ionization densities at such great radii; it is certainly un- 
safe to assume on the basis of present evidence that 
those meteors with large diameter glowing trains may 
not also exhibit ionization columns of comparable size.”* 

The contrary tack was taken by the investigators at 
the University of Manchester. In reviewing the theory 
of meteoric ionization, Herlofson was led to estimate 
that the energy available for the production of ioniza- 
tion would be of the order of 10~‘ of the total available 
kinetic energy.”’ On the basis of these figures, and as- 
suming the ionized molecules were those of atmospheric 
air, he gave a tentative estimate of the relationship be- 
tween visual meteoric magnitude and line density of 
ionization along the trail. A 6th magnitude meteor was 
assumed to produce 10!° electrons per centimeter —of 
path, and a +1 magnitude meteor 10! electrons per 
cm. The electron densities resulting from Herlofson’s 
estimate are so low that an incident wave will not be 
appreciably attenuated by the charges; the received sig- 
nal strength is that caused by Thompson scattering of 
the individual electrons in the trail. A scattering for- 
mula expressing these relations was given by Lovell and 
Clegg who used it to estimate electron densities from 
received signal strengths.?4 They reported a confirma- 
tion of the electron density-magnitude relations pre- 
dicted by Herlofson, but it appears that something was 
wrong with the results. More recent signal strength 
measurements, both of the Manchester and Stanford 
groups have shown that the electron densities are many 
orders of magnitude greater than those given by Lovell 


19 A. M. Skellet, “The ionizing effects of meteors,” Proc. I.R.E., 
vol. 23, pp. 132-149; 1935. 

20 J. A. Pierce, “A note on ionization by meteors,” Phys. Rev., 
vol. 59, pp. 625-626; 1941. 

21C. C. Trowbridge, “Physical nature of meteor trains,” Astro- 
phys. Jour., vol. 26, pp. 95-116; 1907. 

2. W. Allen, Jr., “Reflection of very-high-frequency radio waves 
from meteor ionization,” Proc. I.R.E., vol. 36, pp. 346-353; 1948. 
(See also discussion of this paper by L. A. Manning and O. G. Vil- 
lard, Jr., Proc. I.R.E., vol. 36, pp. 1255-1257; 1948.) 

*8 J. Herlofson, “The theory of meteor ionization,” Phys. Soc. 
Rep. Prog., Phys., vol. 11, pp. 444-454; 1948. 

4 A.C. B. Lovell and J. A. Clegg, “Characteristics of radio echoes 
from meteor trails: I. The intensity of the radio reflections and 
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and his co-workers.18> It now seems that the experi- 
ments of Lovell were in error. In addition it seems that 
the tentative estimates of Herlofson, with which the ex- 
periments were purported to be in agreement, were 
based on rather a weak foundation. They should be re- 
vised using further knowledge of the processes involved. 

What we now know of the structure of meteoric ion- 
ization columns depends on combined inference from 
polarization, reflected signal strength, and duration 
data. Polarization studies show unequal reflection co- 
efficients for parallel and transverse polarization during 
the early part of the echo, as in Fig. 3. No polarization 
effect is ordinarily observed at the end of the echo. Two 
ways of measuring polarization have been used. Clegg 
and Closs transmitted energy with the electric field 
both parallel and perpendicular to the trails of shower 
meteors of known orientations, and measured the two 
reflection coefficients directly.6 Van Valkenberg of the 
Stanford group employed circularly polarized transmit- 
ting and receiving antennas and measured the ellipticity 
of the received echo when the trail was excited with cir- 
cular polarization.’? The method used at Manchester has 
the advantage of identifying the two polarizations with 
the trail axis; the method employed at Stanford enables 
the relative magnitudes of the reflected polarized com- 
ponents to be found for meteors of unknown orienta- 
tion. 

Herlofson was the first to suggest that the cause of the 
polarization effect was a plasma resonance phenome- 
non;”° the resonance takes place when a trail with di- 
electric constant about minus one is excited crosswise. 
The high ‘electron concentration required to support 
plasma resonance is a first clue to the state of the trail. 
A second is that the breadth of the trail must be small 
in wavelengths for plasma resonance to occur; other- 
wise the time lag across the column results in out-of- 
phase excitation of the effect. 

A general theory for the reflection of radio waves from 
meteoric ionization may well start with the assumption 
that the ionization is initially confined within a radial 
region whose scale is determined by the mean free path 
of the air—10 cm or so. If the ionization is assumed 
formed either of air by collision, or of meteoric mole- 
cules, such a distribution would be expected on theo- 
retical grounds. The values of the mean free path are 
consistent with the small size demanded for the polar- 
ization effect. Both Kaiser and Closs,2” and Eshleman!8 
have worked out the reflection coefficient to be expected 
from such a column, assuming an initial Gaussian dis- 
tribution of ionization, and subsequent growth of the 
trail by diffusion. Fig. 6 shows how a Gaussian distribu- 
tion expands due to diffusion alone. Eshleman has con- 
sidered the effect on the column growth of a number of 
values of recombination coefficient, as well as of diffu- 
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DISTANCE FROM TRAIL AXIS — ARBITRARY UNITS OF LENGTH 


Fig. 6—The effect of diffusion upon an initially Gaussian distribution of electrons, When first formed, most 
of the ionization is found within 10 cm or so of the axis. 


sion coefficient. In both cases the calculations were car- 
ried out by finding approximate expressions valid in 
limited size and line density regions. Feinstein?* and 
Herlofson?® had both previously considered the small 
radius low density trail, but they were misled by Lov- 
ell’s erroneous experimental results into believing that 
no need existed for a theory involving greater line densi- 
ties. Herlofson was able to show that in the low density 
region, where no more than 10” electrons exist per cm of 
path, the echo strength would decay exponentially, with 
a time constant determined by the value of the diffu- 
sion coefficient and the square of the wavelength. Such 
a behavior agrees with the observed characteristics of 
many meteor echoes. The mean rate of decay of echo 
strength suggests a diffusion coefficient of about 3 
meters? per second for the heights near 95 kilometers. 
At Stanford for frequencies of 6 to 30 mc, and elsewhere 
for higher frequencies, the dependence of echo duration 
on the square of the wavelength has been checked. 
Not all radio reflections have roughly exponential 
rates of signal strength decay; some show greatly in- 
creased durations and slow loss of signal strength (but 
with violent fading). It has long been realized that the 
great durations could be explained by large meteors 
with great ionization productions, if only recombination 
effects could be neglected. However, it has also long 
been recognized that the recombination coefficient of 


28 J. Feinstein, “The interpretation of radar echoes from meteor 
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10-® cm™ sec! which is appropriate for the H-region 
would quickly reduce even an infinite initial electron 
density to a very low value; the increased durations 
would no longer be possible.*° Computation of the re- 
flections from large highly ionized trails has been carried 
out by assuming that the trails are equivalent to metal- 
lic cylinders; their radius is equal to the radius at which 
the dielectric constant of the actual trail is zero. As has 
been shown by Kaiser and Closs, and by Eshleman, the 
behavior of the echoes computed in this way is quite 
consistent with the observed phenomena of radio re- 
flection if recombination is small, and neglecting for the 
present the fading effects. 

When recombination is included, Eshleman has shown 
in detail the effect of variation of recombination coeffi- 
cient upon line density.!® From his results it follows 
that the maximum coefficient consistent with the longer 
duration echoes is of the order 10—!° cm~ sec7!. Since 
the H-region recombinaton coefficient is 100 times 
greater than this, it may be supposed either that as sug- 
gested by Eshleman, the greatly increased temperature 
along the trail reduces the recombination coefficient of 
air during the initial expansion, or else that in meteoric 
columns it is not mainly the oxygen which is ionized. 
The latter possibility is of especial interest. The result- 
ing recombination coefficient would be determined by 
the meteoric elements, and since we don’t know what 
their coefficients are, we may speculate that they are 
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suitably low, at least for some constituents. Note that 
it is found spectroscopically that the lines observed in 
meteoric light are not those of air; indeed, concentra- 
tions of ionized calcium in particular are known to exist 
in roughly adequate amounts.*? 

Assumption of the ionization of meteoric molecules 
with lower energies of ionization than air may remove 
much of the difficulty in reconciling the kinetic energy 
of the particle with the fraction thereof which can be 
turned to ionization. There are amply enough molecules 
of meteoric material to supply the observed electron 
densities. It is clearly of great importance that experi- 
ments be devised to actually measure recombination 
in meteor trails, and to discover the nature of the ionized 
particles. 


WIND MEASUREMENTS 


It is in the study of upper atmosphere winds that 
meteoric ionization trails are most nearly describable as 
tools used in the process of research. For purposes of 
wind measurement the ionization may be looked at as 
being a marker which makes it possible to follow the 
motions of the air in which it is imbedded. 

The principal means of investigating wind motions 
makes use of a Doppler frequency shift called the “body 
Doppler” ;3? it occurs after the down-whistle and dif- 
fraction fluctuation sketched in Fig. 5. Unlike the vari- 
able frequency and high pitched meteor whistle which 
ushers in the formation of a new trail, the body Doppler 
is of a low and constant frequency for the duration of the 
echo. Five cycles per second is a typical frequency on a 
20 Mc carrier. By scaling the value of the body Doppler 
shift, it is possible to find the radial velocity component 
of the wind drift for a particular trail. It is necessary, 
as Manning, Villard, and Peterson described in the orig- 
inal paper,®? to have a sense indicator which will indi- 
cate whether the motion of the reflector is toward or 
away from the observer, if maximum use is to be made 
of the effect. 

With the aid of the body Doppler shift, it is possible to 
find several quantities descriptive of the wind. The 
speed and direction of the average translation of the air 
mass, including a vertical component if present, is one 
such quantity. This measure of the vector average of 
the individual motions which exist in a region will 
clearly be less than the peak wind speed found at any 
point or in any stratified layer included within the re- 
gion. Another descriptive quantity is the root-mean- 
square wind speed for both horizontal and vertical 
components, as averaged over a height region.*? This 
rms speed is a good measure of the rapidity with which 
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the air is blowing, since air motions in different direc- 
tions do not cancel in producing the final result. 

Average winds are found by a sampling technique; 
the radial velocity component is measured for a large 
number of echoes distributed in azimuth about the ob- 
servatory. In those directions from which the average 
wind is coming, the average radial velocity component 
will be towards the observer, and conversely in the op- 
posite directions. By plotting the average Doppler shift 
versus the azimuthal direction from which the meteors 
are received a sine relationship is found, marred only by 
randon sampling errors. Careful statistical analysis 
shows that the amplitude and phase of the sine varia- 
tion may be used to compute the magnitude and direc- 
tion of the average wind with good accuracy. A sample 
of several hundred meteors must be employed; one or 
two hours of recording generally suffice to obtain a sam- 
ple of this size if sensitive equipment is used in the vi- 
cinity of 20 mc. Velocities found by this means typically 
are of the order of 25 mc, when averaged for several 
consecutive hours over the 80-120 km height range. 
Great variability occurs in both the magnitude and 
direction of the velocity. Average vertical drift is no 
greater than 3 mc as observed at Stanford, and may be 
nonexistent. 

Root-mean-square winds are best found by plotting 
the average square of measured radial velocity com- 
ponents versus the quantity (h/R)*, where h is the 
height of the reflector, and R is its slant range. As has 
been shown by Manning, Peterson, and Villard,?4 when 
plotted in this way the resulting curve is straight, and a 
regression line is easily fitted. The rms horizontal and 
vertical velocities are then found as the (h/R)?=0 and 1 
ordinates. Experimental observations need only be 
taken in the middle range of elevation angles, not over- 
head or on the horizon. Typical horizontal rms velocities 
determined in this way are 50 mc; they exhibit much less 
variability than do the mean vector winds, which appear 
to be averages between large and opposing motions The 
vertical rms velocities are much less than the horizontal 
velocities. They are certainly less than 20 mc on ordi- 
nary occasions, and quite possibly much smaller. The 
resolution of the measurements has not been increased 
enough to give a more accurate answer. 

One of the great advantages of the meteoric method 
for finding velocities is that the motions of discrete 
marked patches of air are examined. On theoretical 
grounds it has been concluded that there is little pos- 
sibility that the ionization does not move with the air 
at meteoric heights. The collision frequencies for both 
electrons and ions are above their gyro-frequencies for 
most of the meteoric region. It is sufficient, however, 
that the ionic collision frequency be well above the ionic 
gyro-frequency in the region in order that the geomag- 
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netic field be unable to affect drift of the ionization wath 
the wind. It may be noted here that in the related 
problem of the effect of the earth’s magnetic field on the 
coefficient of diffusion, again no effect is anticipated. In 
the literature there is a supposed demonstration of such 
a variation of diffusion coefficient,** but close attention 
to the theory shows that the plotted results do not have 
the predicted form. Attempts to verify the effect have 
been made by McKinley, and by the Stanford group 
without success. 

Further development of the Doppler method of wind 
measurement may well take advantage of the definite 
positions of meteoric reflectors in order to delineate the 
variations which exist in wind structure with height. 
The method has been adopted by Huxley’s group in 
Adelaide, Australia; winds are being recorded auto- 
matically, with complete height information accom- 
panying each record. From such work valuable results 
concerning the stratification of wind layers, and the 
variation of tidal motions with height, may be ex- 
pected. 

In addition to those simply described meteoric re- 
flections which exhibit a single Doppler shifted spectral 
component, there are many received echoes which con- 
sist of two or more reflected waves. Beating of the 
separate frequency components causes them to fade in 
amplitude. It has been shown by Manning, Villard, and 
Peterson that analysis of these signals may advantage- 
ously be carried out by use of a “double Doppler” vec- 
tor echo presentation.!° By treating the received signal 
as a vector, spectrum analysis can be achieved by a 
process of double Fourier analysis of two Doppler chan- 
nels in which the beating voltages differ by ninety de- 
grees. In the case of the fading echoes, the fading rate is 
governed by the difference frequency between inde- 
pendent spectral components. To the extent that all 
spectral components are caused by different wind mo- 
tions at separate points along the trail, analysis of fad- 
ing rates can be used to find the difference in wind veloc- 
ity between neighboring reflecting regions.® Investiga- 
tions at Stanford have revealed in addition that in some 
cases it is possible to determine the displacement along 
the trail between the reflecting centers by noting the 
phase difference between the fading envelopes detected 
at separate ground stations. Analysis of fading char- 
acteristics to determine wind gradients suffers from the 
danger that some of the fading may not be caused by 
wind motions at all. Feinstein has shown that for suffi- 
ciently abruptly bounded electron distributions in an ex- 
panding trail, fading may occur as a result of a kind of 
beating of the waves scattered from the near and far 
boundaries of the column.?8 Eshleman has shown that 
such a fading effect can occur without introducing an 
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normalities,” The Phys. Soc., (reports on progress in physics), 
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atmospheric winds at altitudes of 80 to 100 km,” Jour. Atmos. Terr. 
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unduly great increase in the rate of decay of echo ampli- 
tude from its maximum value.!8 In a number of cases 
meteoric echoes are found to fade with two independent 
and well defined fading rates, the lower frequency fading 
envelope limiting the amplitude of the faster fluctua- 
tions. It taxes one’s credulity to believe that all echoes 
of this nature result from four discrete spectral com- 
ponents caused by wind motion, and neatly paired in 
amplitude. Further investigation of the reality of the 
Feinstein type of fading is definitely required. 


THE CONTRIBUTION OF METEORS TO OBLIQUELY 
PROPAGATED HIGH FREQUENCY SIGNALS 


Very little consideration had been given until recently 
to the possible use of meteoric reflections for the trans- 
mission of signals over an oblique path. Even the 
changes in the reflection characteristics exhibited at 
obliquely incidence by individual trails had received 
very little attention. One of the earliest studies of this 
problem was undertaken by Allen,?? who recorded the 
bursts of signal strength received from very-high- 
frequency FM stations over an extended path. He noted 
the existence at times of received signal strengths for 
which the radius of the equivalent metallic reflector 
must have been several hundred meters. More recently, 
the problem has been attacked theoretically by Eshle- 
man and Manning,*? and experimentally by Villard, 
Peterson, Manning, and Eshleman.** The results of both 
theory and experiment indicate that meteors should not 
be neglected in considering possible transmission mech- 
anisms at frequencies above the normal layer criticals. 

The predominant effect of obliquity upon the nature 
of the reflected signal is to increase the echo duration 
by what may be a large ratio as compared with back 
reflection. In the case of back scattering, the factor 
which normally limits the echo duration is the phase 
opposition of energy scattered from electrons situated 
across the breadth of the trail. As the trail expands as a 
result of diffusion, the phase opposition becomes more 
and more pronounced, and the returned signal dies out. 
In the case of oblique reflection, the same process ap- 
plies, but the geometry is now such that a separation of 
sec @ times that necessary in the back scatter case pro- 
duces the same phase cancellation; ¢ is half the forward 
scattering angle. The effect upon phase cancellation is 
as though the wavelength of the incident wave were re- 
duced by the factor cos ¢. Since it has been demon- 
strated experimentally that the echo duration is pro- 


portional to the square of the effective wavelength, it 


follows that the effect of obliquity upon echo duration 
goes as the square of the secant of half the forward 
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scattering angle. This duration increase factor can quite 
easily be as large as ten or twenty times for meteors 
situated near the middle of a path of one or two thous- 
and kilometers. 

In assessing the importance of meteors to obliquely 
propagated signals, it is necessary to consider also the 
total number of meteors in the reflection process. Over 
an oblique path, a large area is available within which 
meteors may fall and produce usable signals. The same 
limitations as in the back scattering case apply to the 
fraction of meteors whose aspect is correct to reflect an 
echo to a given receiver. Although the total fraction of 
meteors effective for oblique transmission is not espe- 
cially different from that for back reflection, the fraction 
effective is a marked function of the part of the region 
between the transmitter and the receiver in which the 
meteor falls. Since the duration increase factor is also a 
function of position, detailed consideration is needed in 
order to demonstrate the total magnitude of the integ- 
rated signal strength effect at the receiver. Also im- 
portant in assessing the quality of meteor propagated 
signals is the probability that there will be at least one 
trail so situated as to propagate a signal at a given time. 
Statistical considerations based upon the above con- 
cepts lead to the conclusion that at frequencies of the 
order of 20 mc, meteoric signals will be present over an 
oblique path a very large fraction of the time. Experi- 
mental study of such paths has shown that useful sig- 
nals are indeed present as a result of meteoric action in 
this frequency range. At the very high frequencies, how- 
ever, it is doubtful that meteoric effects alone are capa- 
ble of supporting a useful transmission. 


THE CONTRIBUTION OF METEORS TO THE AMBIENT 
IONIZATION OF THE E-REGION 


Theoretical study of the contribution of meteors to 
the ambient ionization level of the E-region is hampered 
by lack of accurate information on the correct value of 
recombination coefficient to be applied to meteoric 
ionization, and on the total number of electrons gener- 
ated by meteors in the various magnitude ranges. It 
seems very improbable that the recombination coeffi- 
cient for meteoric ionization could be greater than 
10-8/cm® sec, and if molecules of meteoric matter are 
ionized, the recombination coefficient may be much 
smaller. Assume the largest reasonable value of recom- 
bination coefficient, and make very rough assumptions 
concerning the total number of electrons generated in 
the region; then we may estimate that an ambient level 
corresponding to a critical frequency of the order of 
several hundred kilocycles might result. With what we 
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now know of recombination any such calculation will be 
crude, and has never been carried through with care. 
Nevertheless, meteoric effects clearly are of importance 
to low frequency propagation. In this field, it is impor- 
tant that further studies be made, but they will need 
to be based upon better knowledge of the numerical 
values of the constants involved in the ionization and 
recombination processes. 


CONCLUSIONS 


Since the first serious study of the subject was under- 
taken, the radio detection of meteors has made impor- 
tant contributions to our knowledge in the fields of 
meteoric astronomy, upper atmosphere physics, and 
radio propagation. The application of radio methods to 
find the orbital constants of the meteoric particles has 
been very successful. There is still much which needs 
to be done before an adequate knowledge of the size 
distribution of meteoric particles is available. Study of 
meteoric echoes has yielded dividends in the under- 
standing of the nature of the atmosphere, and of its 
reaction to suddenly applied charge. An important prob- 
lem remains in relation to the recombination coefficient 
applicable to meteoric ionization; is it that for air, or 
that applicable to the meteoric particles? 

Of great value has been the use of meteoric techniques 
in the study of wind motions in the upper atmosphere. 
In this field important and detailed results can be ob- 
tained. With the aid of height finding apparatus, one 
may investigate the structure of the wind in height, 
and determine many statistical properties characteristic 
of its behavior. 

Only recently have meteors been considered exten- 
sively from the point of view of their utility for com- 
munication purposes. It is found that at the high fre- 
quencies meteors can support useful transmissions at 
times when regular layer transmission is absent. How- 
ever, in some situations the meteoric signals are un- 
wanted, and constitute interference. Very little has yet 
been done to investigate the amount of ambient ioniza- 
tion below the E-region which is caused by meteors. 
It is to be hoped that future study will serve to resolve 
the remaining problems which have arisen in meteoric 
research, and to pose many new ones as well. 
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